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ABSTRACT 
The  analysis  of  scattered  electromagnetic  radiation  is  examined  in 
terms  of  Mie  theory  establishing  the  detailed  form  for  the  scattering 
functions.  These  parameters  are  used  to  define  the  criterion  whereby 
the  composition  and  number  densities  of  an  atmospheric  aerosol  population 
may  be  assessed.  Laser  radar,  or  LTDAR,  is  the  proposed  remote  sensing 
device  and  the  governing  system  equations  are  developed.   The  problem  of 
data  inversion  is  surveyed  with  emphasis  on  smoothing  methods,  statistical 
analyses,  and  iterative  techniques*   A  discussion  of  the  numerical 
stability  of  the  solution  is  also  presented.   On  the  basis  of  the  Mie 
model,  inversion  of  the  data  provided  by  the  LTJDAR  probe  is  given  as  a 
rationale  to  effect  an  estimate  of  the  true  particle  number  distribution 
function  in  an  atmospheric  cell. 
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I.   II!TRCDUCTIO'T 

The  advent  of  sophisticated  optical  devices  in  both  the  scientific 
and  defense  communities  has  pointedly  demonstrated  the  need  to  better 
understand  the  primary  propagating  and  interacting  medium  characteristics 
of  the  earth' s  atmosphere.   The  intent  of  this  paper  is  not  to  deal  with 
the  purpose  of  such  devices,  whether  they  are  atmospheric  pollutant  mon- 
itors or  infrared  threat  detectors.  Rather  the  concern  here  is  the  use 
of  these  devices  in  the  real  world  environment  and  the  limitations  im- 
posed on  them  as  a  result  of  atmospheric  optical  signal  degradation. 
This  degradation  is  the  result  of  manifold  scattering  effects  such  as 
Brillouin,  Rayleigh,  Mie,  Fluorescent  and  Raman  scattering  and  includes 
the  atmospheric  constituents'  differential  absorption  of  optical 
radiation.  £<Q 

The  atmospheric  scattering  and  absorption  centers  include  the  air 
molecules,  suspended  particulates  (aerosols)  and  water  droplets  in  fog, 
rain  or  hail.   The  presence  of  aerosols  strongly  determines  the  trans- 
mission characteristics  of  the  atmosphere  and  may  even  dominate  the 
propagation  characteristics.  See  Table  I  and  figure  1  for  a  comparison 
of  particle  sizes  and  the  range  of  number  densities. 

At  the  small  end  of  the  size  spectrum  are  the  air  molecules  them- 

I  1  o         ° 

selves  with  radii  of  10  ^mth   and  a  concentration  of  10   per  cm"".   The 

2      / 

large  end  includes  rain  drops  with  radii  of  10  to  10^/*m  and  a  concen- 

-5      -2      3 
tration  of  10   to  10   per  cm  .   The  "small"  particles,  defined  by  a 

radius  much  less  than  the  wavelength  of  radiation  used  as  the  probe, 
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obey  the  Rayleigh  scattering  model,  exhibiting  a  cross  section  that 

decreases  proportionally  as  the  particle  radius  to  wavelength  ratio  to 

the  fourth  power.  It  is  this  fourth  power  law  for  molecular  scattering 

that  Rayleigh  used  to  explain  the  blue  color  of  the  daytime  sky. 

TABLE  I 

Particle  tyoes        Particle  Radius  (xm)        Concentration  (cm  y) 

Aitken  nuclei  »         1Cf3  to  10"2  102  to  10^ 

haze  particles  10~2  to  1  10  to  103 

fog  droplets  1  to  10  10  to  102 

cloud  droplets  1  to  10  10  to  300 

-::'  Initial  nuclei  present  in  the  atmosphere  from  the  earth's  surface 
[(J   and/or  meteoric  sources. 

For  wavelengths  larger  than  those  of  visible  light,  this  scattering  is 
very  small.  As  the  particle  size  approaches  that  of  the  radiation  wave- 
length, the  complete  Mie  theory  is  required  to  explain  the  scattering 
characteristics.  The  dependence  on  wavelength  is  highly  oscillatory  and 
very  complex.  As  the  particle  grows  still  larger,  its  behavior  in  the 
radiation  field  approximates  that  of  a  large  solid  object  casting  a 
shadow  field  in  proportion  to  its  geometrical  cross-section.  The  Mie 
theory  includes  both  Rayleigh  scattering  for  small  radii  and  the  "non- 
selective" scattering  in  the  limit  of  large  particle  radii.  The  theory 
is  covered  in  some  detail  in  the  next  section  of  this  paper. 

Once  one  has  both  the  atmospheric  density  and  data  on  the  aerosol 
composition  and  size  distribution  as  functions  of  altitude,  location, 
and  time,  the  transmittance  can  be  computed  on  a  real-time  basis  by  some 
rather  sophisticated  computer  models. 

Aerosols  are  colloidal  particles  dispersed  and  suspended  or  falling 


10 


slowly  in  the  gaseous  mixture  of  the  atmosphere.  Some  examples  are 
smokes,  haze,  clouds,  small  droplet  (less  than  1/cm)  fogs,  and  finely- 
divided  soils.   These  aerosols  have  as  many  sources  as  they  have  diff- 
erent forms  and  compositions.   The  chemical  reaction  of  various  atmo- 
spheric gasses  whose  minor  components  are  nitrous  oxides,  terpenes,  and 
hydrocarbons  can  produce  solid  particles.   The  presence  of  combustion 
products  from  inefficient  burning,  which  may  be  locally  critical,  and, 
on  a  much  larger  scale,  from  forest  fires  and  volcanic  action  serves 
only  to  increase  the  atmospheric  load  of  particles  that  may  take  years 
to  settle.  Over  large  bodies  of  water,  the  dispersion  of  solutions  into 
the  near  atmosphere  region  brings  about  the  formation  of  fogs  and  rain 
as  well  as  clouds.   As  a  class,  aerosols  may  either  help  or  hinder  the 
researcher.   These  particles  may  be  used  as  tracers  enhancing  a  received 
signal  to  detect  atmospheric  motions,  distributions,  and  chemical  com- 
position.  However,  they  primarily  provide  a  large,  random  source  of 
environmental  noise.   In  this  latter  mode,  the  particle's  presence  tends 
to  confuse  the  received  signal  at  the  sensor  and  provide  noise  that 
originates  outside  the  system  electronics,  making  it  difficult,  if  not 
impossible,  to  supress.  jjoj  Historically,  aerosols  have  been  divided 
into  classes  to  ■  facilitate  their  theoretical  modeling.  These  models 
involve  several  free  parameters  that  are  designed  to  "best  fit"  the 
observed  data.   The  classes  are: 

1.   Maritime  -  This  group  is  characterized  primarily  of  sea  salt 
with  some  water  of  hydration  and  a  chlorine  concentration  slightly  higher 
than  sea  water.  Salt  particle  concentration  decreases  rapidly  above  a 
height  of  about  500  m. 
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2.  Rural  -  The  largest  number  of  components  are  silicon,  iron 
and  some  sulfates  with  approximately  10%  to  15%   organics.  Of  this  total, 
about  l/3   of  the  aerosol  is  hygroscopic,  whereby  air  moisture  is  readily 
absorbed.   This  same  class  is  also  found  at  sea  at  higher  altitudes  with 
a  large  particle  density  falling  off  rapidly  as  distance  at  sea 
increases. 

3.  Urban  -  Combustion  products  and  the  issue  from  industrial 
processes  make  up  the  bull:  of  this  class. 

A.   Troposphere  -  This  is  an  atmosphere  division  lying  above  the 
surface  boundary  layer  and  consists  of  the  rural/maritime  classes  with- 
out large  particles. 

5.  Stratosphere  -  Primarily  (30^  to  90%)  consisting  of  the  large 
sulfate  ions  and  particulates,  this  group  also  shows  a  periodic  large 
influx  during  volcanic  activity. 

The  typical  atmospheric  analysis  technique  is  to  use  two-ended  sys- 
tems requiring  a  source  of  radiation  of  known  properties,  a  propagating 
path  through  the  atmosphere  of  specific  length  and  reasonably  known 
conditions,  and  a  stable  receiver  at  the  path's  end.  Experiments  such 
as  these  are  able  to  measure  aerosol  size  and  composition  distributions 
to  check  out  the  various  composition  and  scattering  models.   They  also 
predict  the  performance  characteristics  of  a  particular  optical  device. 
However,  two  well- separated  platforms  are  required  with  the  attending 
problems  of  phasing  and  communications.   Also,  only  low  altitude,  hori- 
zontal work  can  be  done.   A  one-ended  system  would  eliminate  the  need 
for  one  of  the  platforms  and  enable  a  full  altitude/azimuthal  capability 
in  a  real-time  environment.  There  is  a  system  that  can  easily  be 
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adapted  to  the  one-ended  process  and  is  available  today.   The  laser 
radar,  LIDAR,  system  is  similar  in  operation  to  Radar  with  optical 
wavelengths  instead  of  radio  waves.   In  fact,  the  acronym  is  somewhat 
the  same,  meaning  "Light  Detection  and  Ranging".  Conceptually,  a  laser 
pulse  is  emitted  into  the  medium  to  be  analyzed.  This  pulse  interacts 
with  the  components  of  the  medium  and  is  both  scattered  and  absorbed  in 
distinct  ways  depending  on  the  material  type  and  density.  Some  portion 
of  this  pulse  is  returned  as  an  echo  to  the  point  of  origin  where  a 
telescope-like  receiver  is  coupled  to  a  photcdetector  that  converts  the 
received  optical  radiation  into  a  proportional  electrical  signal.  The 
pulse  return  timing  is  gated  to  the  range  of  the  atmospheric  cell  under- 
going analysis.  The  intensity  of  the  echo,  the  echo's  polarisation 
change,  and  any  measured  frequency  shifts  all  undergo  involved  processing. 
The  outputs  can  consist  of  composition,  type,  and  size  distributions. 
Some  state  of  the  art  applications  of  the  technique  include:  metero- 
logical  LIDAR  investigations;  middle  atmospheric  LIDAR  studies;  tropo- 
spheric  chemistry  and  diffusion  research;  studies  of  atmospheric 
propagation  and  radiative  transfer;  absorption,  Raman  and  fluo.rescent 
spectroscopic  applications;  and  pressure-temperature  measurements  of 
atmospheric  layers.  £2}  A  recent  theoretical  paper  indicated  that  a 
series  of  measurements  involving  optical  scattering  and  absorption  can 
yield  about  7  or  8  independent  pieces  of  information  concerning  the 
aerosol  size  distribution  within  the  analyzed  cell,  covering  a  range  of 
particles  with  radii  from  O.l^mto  5/<ra.  £2CJ 

Consider  such  a  single-ended  LIDAR  analysis  system  on  board  the  type 
of  mobile  platform  in  use  today  by  both  scientific  researchers  and  the 
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military.      The  optimum  performance  of  a  specific  optical   scanning/de- 
tection/tracking device   could  be  evaluated  in  a  real-time  scenario   with- 
out the  need   to  resort  to   generalized  models   applicable  only  in 
standardized  areas. 


u 


II.   GENERAL  SCATTERING  THEORY 

Scattering  is  the  process  whereby  an  incident  electromagnetic  wave 
loses  energy  and  some  portion  of  this  energy  undergoes  re-emission.   In 
general,  the  scattered  wavefront  has  a  different  propagation  direction 
than  the  incident  wave.  See  figure  2. 

For  particles  with  a  low  number  of  electrons,  the  scattering  can  be 
computed  using  quantum  mechanics.   In  this  sense,  the  incident  wave  can 
only  interact  with  the  stationary  states  of  the  scattering  particle's 
electronic  system  within  a  small  bandwidth  of  frequencies  about  the 
allowed  transition  levels.  One  can  discuss  Rayleigh  scattering  as  the 
absorption  of  the  correct  energy  photon  and  the  subsequent  re-emission 
of  a  photon  of  essentially  the  same  energy.  An  energy  level  diagram 
would  look  like  figure  3. 

Likewise,  Raman  scattering  involves  a  definite  frequency  shift  be- 
tween the  incident  and  scattered  radiation.   If  the  scattered  frequency 
is  less  than  that  of  the  incident,  the  process  produces  the  Stokes 
lines  spectroscopically.  And  when  the  reverse  is  true,  the  anti-Stokes 
condition  exists.  C&3.  See  figure  4-» 

However,  for  molecular  and  larger  particle  scattering,  the  collection 
of  electronic  states  becomes  great  enough  to  let  us  deal  with  the  prob- 
lem classically.  The  particle  is  considered  to  be  a  solid  dielectric 
with  a  complex  index  of  refraction. 

One  of  the  most  widely  studied  problems  in  diffraction  theory  is 
electromagnetic  and  acoustic  wave  scattering  by  a  sphere3   The  problem 
of  diffraction  of  a  plane-monochromatic  wave  by  a  homogeneous  sphere 
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characterized  by  its  radius  and  complex  index  of  refraction  was  rigor- 
ously solved  by  G.   Mie  in  1908.      His  analysis  involved  the  solving  of 
Maxwell's  equations  for  the  electromagnetic  field  when  an  incident  plane 
wave  interacts  with  a  material  interface.     V/hen  the  discrete  boundary 
separates  media  of  different  optical  characteristics,    a  scattered  wave 
is  generated. 

Today,    an  extensive  bibliography  of  approaches  to  this  problem  exists. 
The  author  will  follow  the  analysis  presented  by  Kerker  (1969)    Q 2^ 
throughout  this  section  in  order  to  present  a  consistent  derivation. 
First,    the  salient  features  of  the  Maxwell   theory  will  be  listed,   lead- 
ing to  the  wave  equation.     The  solution  satisfying  the  appropriate 
boundary  conditions  will  be  obtained;    the  resulting  infinite   series 
solution  not  only  completely  defines  the   scattered  wave,    but  gives   the 
electromagnetic  states  of  the  particules1    interior  as  well.      The  single 
particle  scattering  results  will  be  compiled  and  extended  to   cover 
collective  scattering. 

Maxwell's   equations,    in  rationalized  MKS  units:    Q2j,0ij 
Coulomb's  law  ^  •  Q      =    (°  (1) 

Ampere'  s  law  p  x  jj      =  —  fc  t  J  (2) 

Faraday's  law  V*£     -Tf^  (3) 

absence  of  free  magnetic  poles  V  «B  -q  (4-) 

The  symbol  J£  denotes   a  vector  quantity.      Following  are  definitions 
of  vectors  used  in  derivation: 

E  -  electric  field  intensity  (volts/meter) 

H  -  magnetic  field  intensity   (amps/meter) 

o 
D  -   electric  flux  density   (dielectric  displacement,    coulomb/meter"" 


is 


o 
J3,  -  magnetic  flux  density  (magnetic  induction,  weber/meter  ) 

2 
-7  -  electric  current  density  (amps/meter  } 

3 
/°  -  electric  charge  density  (coulomb/meter  ) 

t  -  tine  (seconds) 

These  equations  will  define  the  electromagnetic  field  configuration 

for  all  points  in  space  and  uniquely  determine  the  field  vectors  giver. 

the  media  characteristics  below. 

fe -""!!,  (6) 

J-'^  (7) 

The  media  parameters  thus  defined  are: 

&    -  permittivity  (electric  inductive  capacity,   farads/meter) 
yu<   -  permeability   (magnetic  inductive   capacity,    henrys/meter) 
&    -   conductivity   (specific  conductance,    1 /ohm-meter) 
In  free  space    (vacuo): 

£*  £  =  S. 35-42  X   10~12  farad/m 
^s^-  l&l  10"     henry/m 
As   a  consequence  of  the  invariance  of  the  Maxwell   equations,    the  quantity 
{f&T~   has   the  units  of  velocity  and,    again,    in  free  space: 

C*   (>.$.)"*    =   2.W125    X     io8    „/sec  (8) 

This  last  quantity  is  the  speed  of  light  in  vacuum.      When  the  quantities 
€,j*t    and  o^ are  independent  of  direction,    the  region  in  space  they 
characterize  is   said  to   be  isotropic.      Two  more  equations   ore  required 
to  complete  the  picture, 

Lorentz  Force  £  =  p  (H  +  U   X  §)  (9) 

Newton's   second  law  of  motion     £  ~     mo.  (10) 
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'■/here:  ?  -  force  on  charged  oarticle  (nevtons) 

u  -  velocity  of  charged  particle  (m/s) 

m  -  mass  of  charged  particle  (kg) 

a  -  acceleration  of  charged  particle  (m/s  ) 
All  macroscopic,  linear,  electromagnetic  phenomenon  are  described  by 
equations  (1)  through  (10)  with  the  associated  quantity  definitions. 
±'ne   classical  dynamics  of  interacting  charged  particles  and  electro- 
magnetic fields  are  thus  specified. 

A  certain  degree  of  convenience  is  obtained  using  the  definitions  of 
the  two  following  auxiliary  vectors: 

£  =  £-*.£  (ii) 

m  =  0£Ojs-h  (12) 

Where:  ?  -  electric  polarization  (c/m  ) 

M  -  magnetic  polarization  (A/m) 
One  can  see  that,  in  free  space,  these  quantities  vanish  identically  and 
they  may  then  be  considered  as  a  measure  of  the  effect  of  matter  on  the 
local  fields.  Some  dimensionless  quantities  may  also  be  defined: 

K£  =  f£.  (13) 

Km--^//<o  (u) 

'.-.here:     jff    -   specific  Z>  (relative  dielectric  constant) 

W     -  specific/*  (relative  magnetic  permeability) 
The  media  may  further  be  described  with: 

M^ma  (i&) 

Where:  -^e~  electric  su s ceo t ability 

A~-  magnetic  susceptability 
The  previous  set  of  equations  can  be  combined  to  obtain: 
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xm=  yr^-i  (is) 

Maxwell's  equations  apply  to  regions  where  the  media  may  be  repre- 
sented by  the  continuous  functions  C,/-*,  and  CT.  When  a  medium  discon- 
tinuity is  encountered,  such  as  the  conditions  across  a  regional 
interface,  the  fields  must  satisfy  a  specific  set  of  relations  at  the 

boundary.  See  figure  5.  Let: 

2 
£$   -  surface  current  density  on  interface  S  (A/m  ) 

/f  -   surface  charge  density  on  interface  3  (c/m  ) 

The  Maxwell  equations  (1)  -  (4-)  are  cast  in  point  form  for  each 

differential  point  in  space.  One  can  apply  various  integral  theorems 

to  change  the  same  equations  to  integral  form  and  use  them  to  define  the 

conditions  that  the  fields  must  satisfy  by  crossing  the  interface.  JJ7J 

The  boundary  conditions  obtained  are  inherent  in  the  Maxwell  equations. 

Boundary  conditions : 

(a)  the  tangential  component  of  E  is  continuous 

(£2-£»)x  £  =  O  (19) 

(b)  The  tangential  component  of  H  undergoes  a  discontinuity  eaual  to  J 

f*  x       <~s 

(M»-Hi)*  a=  Is  (20) 

(c)  The  normal  component  of  D,  undergoes  a  discontinuity  equal  to  A 

(d)  the  normal  component  of  B  is  continuous 

(£*-§.)•  a  =  o  (22) 

For  all  cases  of  dielectric  scattering,  it  is  assumed  that  there  is  no 
free  surface  charge  or  current  density.  This  assumption  begins  to  break 
down  for  large  particles  in  the  presence  of  storm  activity.  However,  we 
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will,  for  simplicity,  assume  that  J  -  0  and/0  =  0.  Therefore,  the 
tangential  components  of  S  and  H  and  the  normal  components  of  D  and  B 
are  continuous  across  the  boundary. 

Charge  conservation  is  included  in  the  Maxwell  theory  and  is  written 
as: 

continuity  equation   V  '  5v  +  "tfN  =  ^  (23) 

Where:  J  -  volume  current  density  (A/m   ) 

■a 

xO„  -  volume  charge  density  (c/m  ) 

Throughout  the  following  derivation,  a  few  vector  identities  will  be 
required. 

Ilote:  A  -  continuously  differentiable  vector  field 
0  -  continuously  differentiable  scalar  field 

V  •  (V  X  A  )  =  O 

r  x  (^  0)  =  o 

v  x  (?xd)  =  vCv-4")  -  V2A 

if  the  fields  described  by  (1)  -  (4)  are  invariant  under  a  gauge  trans- 
formation, then  they  are  well  defined  by  the  preceding  formalism. jj6J 
We  define: 

A  -  vector  potential  field 

0  -   scalar  potential  field 
in  the  following  relations. 

B  -  V  XA  (25) 

Compare  this  to  equations   (4)   and  (24.).     Using  equation  (25)   in  (3),   we 
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(24) 


get  the  following: 

How,   refering  to   (24), 

=>    I+|c4-    -*rf  .thus 

A  and  0  are  not  unique.     Let  g  be  any  scalar  function,   then  the  new 

potential  can  be  defined  as: 

&'   =     A  -  V9  (27) 

jZT.^  +  ^9  (28) 

Now  the  fields  E  and  3,  are  recomputed. 

B'  =  r  x  (d-^9)  =  B  -tx  (79}  '-  5 
=>    B'  =  8 

=>    £'=£ 

The  physical  observables  are  the  fields  E   and  3  and  it  can  be  seen 
that  the  gauge  transformation  defined  by  (27)  and  (28)  leaves  these  fields 
unchanged.  Thus,  the  Maxwell  equations  are  gauge  invariant. 

We  now  turn  from  formalism  to  a  more  practical  aspect  of  electro- 
magnetic theory  and  introduce  the  Wave  Equation: 
PX(i):  n(r)(E):  *|(tXB) 
Note  that  the  time  and  space  operations  of  differentiation  are  commutative 
for  these  continuously  differentisble  functions. 

00,(24):  v  (▼•£)-  v\  =  -^r^x^ti/) 

The  assumption  is  always  made  that  all  space  and  time  derivations  of  the 
parameters  6,aa,  and 0  are  either  slowly  varying  with  respect  to  the  rate 
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of  change  of  the  fields  themselves  or  identically  equal  to  zero. 

CO,  U),  (5),  (7)  :  v-  (l  £  )  =  £  V  •  E  =  P 

The  same  argument    also  applies  to  the  quantity  i°/ 6  concerning  its  rate 
of  change  with  the  rate  of  change  of  the  fields. 
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St    ~    '       '     3t: 


Therefore:   V  *  fe  -  ^^^  £   '  ^"|f^  &  -  °  ^ 

This   equation,    (29),    is   the  well-known  nonhomogeneous   equation  for  damped 
wave  motion  -   the  ',/ave  Equation.     Similarly,  V  X   (2): 

vx  Cvx  t!)  =|t(v^5^Crxj) 

CO,  (0,(2*0:  v(T7.i4Vv1a=-^lr^xS)+^^xE) 

CM)    -=>  V  •  (•«■  K  0  -  >"  V-  (i  =    0=r>    V  •  K    =  O 

->  *V  <"*&*- £Hf»*=0  (30) 

The  relations  for  non-dissipative  media  are  obtained  by  setting  &  -   0. 
Both  equations  (29)  and  (30)  actually  represent  3  equations  each.  Each 
component  of  the  vector  field  must  satisfy  these  equations  for  the  vector 
field  itself  to  do  so.  Let  the  symbol  u"  correspond  to  any  of  the  six 
components  of  the  two  fields  E   and  H.   The  scalar  wave  equation  can  be 
expressed  by: 
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The  fields  described  by  (31)  represent  a  transfer  of  energy  in  both  space 
and  tine.     The  electromagnetic  wave  has  a  time  periodicity  that  sug- 
gests a  form  for  the  energy  flux  crossing  a  unit  area  per  unit  time: 
S    =  E  x  H  (32) 

Equation  (32)   is  the  Poynting  vector  equation,   whose  time  averaged 
magnitude  over  the  period  of  the  wave  yields  the  wave  intensity, 

I  =    VVJj  Si  dt  (33) 

where T is  the  period  of  the  wave.  The  field  symbols  E  and  H  in  (32)  are 
the  real  parts  of  the  complex  fields  whose  components  satisfy  (31). 

The  form  of  the  scalar  wave  equation  suggests  a  harmonic  time  de- 
pendence of  the  form: 

where  o>  is  the  angular  frequency,  i  =  V(-1 ),  and  r  is  the  usual  three- 
dimensional  position  vector.  Note  that: 

au      •    n  a^a       .   2.,, 

yf    -  lcoU.  and     gTT"     =     ^    ^ 

then,    (31):      72U  -  ^£a;  ^  d^a;^    =0 

=>   V2U  +  (/**  £,  to1  -  [/a  cr  co)  U  =  O 

Define  the  propagation  constant  k  by: 

K2  =  /^t^1-  i^auj  (35) 

then:  v    U.  *■    K2  U    =    O  ^36^ 

Equation  (36)   is  the  Helmholtz  Wave  Equation,  Another  form  of  (35)   is: 

K  =    <xx   -  c/3,  (37) 

where:      *,  B    -ftStf/*^  V  /J}  *  (») 

*■  =  4f/r/*72-'F 

It  is  the/3.  term  that  determines  the  damping  undergone  by  the  wave  per 
unit  path  length;  the  energy  is  dissipated  as  Joule  heat  in  the  medium 
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whentf  >  0. 

Recalling  that  the  quantity 

^  =  C^)*  (40) 

is  a  velocity  which,  in  free  space  is  c,  the  relation  for  an  electro- 
magnetic wave  in  a  medium  described  by  light  speed  V  is 

where:  \  -   wavelength 

v  -  frequency 

V    -  velocity  of  the  wave 
Also,  ZVf-0  r  a)  (42) 

and  this  is  the  angular  frequency,  uJ ,   mentioned  earlier.      In  free  space, 

The  index  of  refraction,   n,   is  defined  as 

*  =  c/v  (u) 

(41),    (43),    (M):        r»a    ^J-  =       %. 

=>       n  A   =  XQ  (45) 

It  is  clear  that  A6  is  the  wavelength  when  the  wave  is  propagated  in 
free  space.   In  the  same  medium,  Cf  -   0  implies  that 
/6,  a  0   =>"f,a  u)  (/*«  £a)  2  =  u>/c 

=>      K0  .  lo/c  z:    2lTv>A0V 

=>      K0=    2TT/^o  (46) 

This  k     is  the  free  space  propagation  constant. 

The   complex  refractive  index,   m,    can  be  defined  as 

m  =  k/k0=  n  0  -iVO  (47) 

where:  ">C  -  index  of  absorption  (extinction) 

As  an  example,  suppose  that  there  is  a  plane  wave  propagating  along  the 
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z  -   axis  whose  x  -   component  is": 

\e 


■=>fcx  =  fle  e 

— nl"  X  z 
The  attenuation  of  the  wave  is  seen  explicitly  as  e  "  "°  .  An  expansion 

of  this  example  shows  that  the  resultant  plane  wave  solution  leads 

directly  to  the  Stokes  parameters.   Assume  the  propagation  direction  is 

z,  as  stated  earlier,  and  that  the  monochromatic  wave  is  traveling 

unbound  in  a  source  free  isotropic  medium.   It  is  assumed  that  the  wave 

car.  by  completely  specified  in  the  xy-plane  b-j   the  z,t  coordinates  of 

the  wave.     Thus: 

v2u  =  iluCr.n    =s>  au    a* 

(34),   (36):  gp^-  f  «lto]eLfll*a0 

^-  f      ,.  2  0 

bz 
Substituting  a  trial  solution,  f  =  ae  , 

=>  a.l>V%  K1aaki=  o       =>    br=  -k%  t-*tk 

From  the  definition  of  k,  -ikz  is  selected  as  the  value. 

This  equation  describes  the  form  for  each  of  the  six  vector  components 

of  E  and  H.   Recalling  (1),  under  the  oresent  constraints: 
3£y     <^cu     ^£7 

(1)     T7"  *    T4     +■    9l      =  O  and    a=    £x?^  )ancI  ££. 

But,   ^a      ^a     _   _s  I^zl      ^  ^.^   c     -  n 
*=>£,«  4«'***"*x>     a**   £^=  *ye 
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(3!' (6)    vxf=  -<H±a  _>      -<aULx  .   --^y   , 

at  "  <j*  *  -     9z 7    ' 

it  is  quite  sufficient  to  describe  the  behavior  cf  3  as  the  value  of  II 

can  be  readily  obtained  using  equations  (l)  through  (4.).  Given  the  form 

of  u(s,t): 

j^-X  =  ia>Ht        ,         yf  Y    =    LU}tfy 

-K,  _    r 
•=•>      "V*«  £  CO  My     =.     ~  c/C  Cx  —S'       H  kj     s    yu  U>         * 

Note:    £x/h(y    =     ~EHjH%    =/><u>/k 

The  terrn/^d/k  can  be  reduced  using   (35)    to   the  following: 

Defining  the  intrinsic  imoedance  of  the  medium  for  Diane  waves  as 


Zo    -     (<5-^J    l  US) 


u> 


Finally,  £.for  the  plane  wave  is: 
t(i0*  -Kz  ^^ 

Ex  =   Ae 

ez=  0 

where  <J    and  d^are  arbitrary  phase  angles   to   account  for  the  general  his- 
tory of  the  wave.     The  II  wave  is  also  determined  by  equation   (4-9): 
I  i(u>t-K2.  v-oT2) 

H  ,«£,,.  <&*-«**'■>  r  (50, 

Hz  =    O 

l'/e  see  that:  _    ,,        _    ,,  -e   c  PC 
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=>         E«H  rO  (51) 

Thus,  3  and  H  are  perpendicular  to  each  other  for  the  plane  wave 

solution. 

The  concept  of  polarised  light  can  now  be  introduced  with  the  help 

of  equation  (49),$,,  and<%.   The  general  plane  wave  solution  describes 

an  elliptical  rotation  of  the  1  vector  in  the  olane  oeroendicular  to  the 

direction  of  propagation.   Defining  the  phase  difference  asJ": 

S   a  Sz  -  <f,  (52) 

The  representation  of  the  state  of  a  polarized  wave  can  then  be  uniquely 

determined  using  the  Stokes  parameters,  which  are: 

5e=   A2+82  =   flf  +   B?  (53) 

S,=  42-E>2  =  S0cos2  Vcos  2  X  r5/} 

S2=    2  /)3  cos  /=  S0SJrtZ  V  Cos  a  X  (55) 

S3-   2#6  SlNcTr  S0S|N2  X  (56) 

The  lengths  A-    and  3,    and  the  angles  r,    and X  are  illustrated  in  figure  6. 

The  quantity   (u>t  -  kz)    can  be  eliminated  from   (49)   using  equation 

(52),    where  E     and  E^r  refer  to   the  real  parts  of  the   complex  wave: 

(f£Y  IfeY-    2ix£Y  COS  is    SIN2  S  (57) 

VAX       \  B/  A3 

This  is  the  equation  of  an  ellipse.   The  ellipse  in  figure  6  is  dra'.ni 
"right-handed",  that  is,  the  head  of  the  electric  field  vector  follows 
the  arrows  when  viewing  the  wavefront  face  on.   In  order  to  exhibit  right- 
handed  polarization,  sinS>0  and  OiX£ff/U*     Linear  polarization  exists 
when: 

S»jtt,     jso.st,  i2,£3, .-  (53) 

and  circular  polarization  corresponds  to 

S-J&jir,  j  =  tijt3,t5,  ±7,  ...  (59) 
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Pictorial  definition  of  the 
quantities  A, ,   B.  ,  H*  fX 
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The  Stokes  parameters  satisfy: 

S0    -      5  I    +  S2    +   53  (60) 

and  a  complete  state  of  polarization  is  obtained  from  only  three  in- 
dependent parameters. 

Natural  light  is  defined  as  the  state  where  the  intensity  of  light 
is  the  same  in  any  direction  perpendicular  to  the  propagation  vector,   k. 
The  phases  of  the  various  components  in  natural  light  are  randomly 
varying  so  that  no   correlation  exists  between  them. 

=r>        S?    =    Snz    =   S3    =  O  (61) 

The  combination  of  eUiptically  polarized  beams  can  be  accomplished  using: 

%'-£  Vc  ,    S,=  i  S,i ,    S2--  i  52i  ,    S3  --  t  s3l  (62) 

Then  for  non-identical  polarization  states  for  all  beams: 

S0     >     S,     +    Sz    m   S3  (63) 

Partially  polarized  light  can  be  described  as  a  superposition  of  natural 


(64) 


light  and  a  set  of  incoherent  beams: 

Finally,  the  fraction  of  polarization,  F,  can  be  defined  as  the  intensity 
ratio  of  the  polarized  part  of  the  wave  to  the  total  wave  superposition: 

F  »l£r  --  »?/(»!♦  s.f)  Co<f«i)  (65) 

It  can  be  shown  that  the  parameters  of  a  scattered  wave  are  linearly 
dependent  on  the  parameters  of  the  incident  wave.jjoj        This  relation 
is  written  as  follows: 

I,=   (l/R*JgIi  (66) 

where:  R  -  observation  distance  to  scatterer 
i  -  incident  wave 
s  -  scattered  wave 
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I  = 


The  modified  Stokes  parameters  are  defined  first  by  writing  equation 

(4-9)   as:  l^t  'LU)t 

£X=E'e  '      £Y=E2«  (67) 


tnen, 


I.    =    |£,|2  (63) 


ii  =  i£.r  (69) 

U  =     2    Re    iMz*$  ^7°) 

V  ^    2   Im   {Ma*J  (70 

The  4^A  Stokes  matrix  a,  looks  like  this: 


C7  = 


(72) 


where  the  f . .  are  defined  in  terms  of  the  scattering  functions  given  in 
j~24J.  For  spherical  particles,  the  quantities  ^10=  ^21  "  ^  an<^  5^ 
contains  four  constants. 

Having  laid  the  foundation  for  the  Mie  treatment  of  scattering 
phenomenon,  let  us  now  consider  a  plane  electromagnetic  wave  incident  on 
an  isotropic,  homogeneous  sphere  of  arbitrary  size.  The  object  interferes 
with  the  free  propagation  of  the  wave  and  secondary  waves  are  generated 
that  are  predictable  by  Maxwell's  equations  and  the  boundary  conditions; 
and  these  secondary  waves  form  a  scattered  radiation  field.  The  incident 
wave  should  be  linearly  polarized  for  two  reasons.   First,  if  the  light 
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is  unpolarized,  then  it  can  be  resolved  into  two  linearly  polarized 
components  each  of  which  act  independently  of  the  other.   Gne  only  needs 
to  consider  the  interaction  of  the  sphere  with  a  plane  (linear)  polarised 
wave  and  then  apply  the  principle  of  superposition  to  obtain  the  total 
effect.  The  LIDAR  laser  pulse  is  plane  polarized,  by  definition,  and 
this  is  the  wave-type  that  is  considered  for  tropospheric  probing.  A 
geometric  diagram  of  the  problem  is  shown  in  figure  7,  suggesting  a 
spherical  coordinate  system  as  the  most  convenient  for  analysis. 

The  fields  of  electromagnetic  energy  in  the  region  of  space  filled 
by  the  particle  and  the  incident  wave  are  resolved  into  three  components: 

a.  the  incident  wave  -  E.,  H. 

b.  the  Darticle  interior  wave  -  2  ,  H 

T*  +*r 

c.  the  scattered  wave  -  E  ,  H 

*  s  -*s 

For  these  quantities  to  properly  define  the  fields,  they  must  each 
satisfy  Maxwell's  equations  (1)  -  (4-),  the  wave  equations  (29),  (30), 
and  the  boundary  conditions  (19)  -  (22).  The  internal  field  E  ,  H 
must  match  the  external  field  2.  +  E  .  H.  +  H  in  accordance  with 

~1    ^S    ~1    *'S 

these  relations. 

KerkerQ2]et  al.,  suggests  that  the  derivation  deal  with  the  scalar 
wave  equation,  (31).  Two  new  functions  may  be  introduced,  the  electric 
Hertz  vector,]!^,  and  the  magnetic  Hertz  vector,!^,  and  are  defined  as 
follows: 

B,  =y"£tr  x jt-  jTi  (73) 

E,=  v(v-%)-/«ej&  Ui  (74) 

&2=->«£V  x^-TT2  (75) 

M2=  v(v-Tri)-s«Z^  &-<*;<<%.  (76) 
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Scattering  of  a  plane  wave  by  a  sphere 
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Suitable  transformations  may  be  applied  to  define  other  Hertz  vectors 
as  long  as  the  fields  are  left  invariant.  The  vectors  thus  defined  must 
satisfy  the  vector  wave  equation: 

* '3 -*/"£•  !K -£./*&  IT  =  -tZ  (77) 

*  *r.  -  **..&  %-  w&  tr»  =  -  a  (vs) 

where  the£  and  J£  result  froa  distributions  of  electric  and  magnetic 
dipoles.  Therefore,  two  vectors  can  be  used  to  describe  all  the  fields 
present. 

The  transverse  magnetic  wave  (TM),  also  called  the  electric  wave, 
has  a  zero  value  for  the  radial  component  of  its  magnetic  field  intensity, 
H.  =  0.  Likewise,  the  transverse  electric  wave  (TS),  the  magnetic  wave, 
has  a  zero  radial  component  of  electric  field  intensity,  E„  =0.  The 
TM  wave  is  conceived  as  being  the  result  of  oscillating  electric  di- 
poles in  a  spatial  region,  and  the  TS  wave  as  the  result  of  oscillating 
magnetic  dipoles.  The  total  solution  is  a  construct  of  these  two 
wavefronts. 

The  Kertz-Debye  potentials  are: 

"V"2i-^  (79) 

-V  'JTZ  -  TT2  (so) 

These  potentials  are  the  solutions  of  the  scalar  wave  equation.  Solving 
(79)  and  (SO)  for  the  potentials  and  applying  the  appropriate  boundary 
conditions,  one  can  then  find  the  field  vectors  describing  the  TM  and 
TS  waves.  The  total  field  may  be  derived  by  the  direct  addition  of  the 
component  waves. 

The  following  equations  for  the  total  field  vector  components  can  be 
obtained  by  applying  these  concepts  to  the  above  relations.   In  spherical 
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coordinates,    we  have: 
£r-Elr+£2r    =  £l   Crir,)+    K  *  C  r  ^  )  +  O  (31) 

He^ie*E29  =7 l7F©^rTr'^K2'Fsm9"a5r  (r>0  (32) 

E^  =  E^  +  £2^  =  FT^i"^a?(frrr1>)-K^-|^C''^  (33) 

Mr  ^   tf»r  +H2r  =  O  +^  (riri)+  k'CrlTj  (34.) 

He=H19  +  H2a^  -K.Tfesi^Cr^^^C'rfrO  (35) 

W^.W,gr  +  H2i/=    K.-jr-^CrTT^  +  Ts^aTa?^^  (36) 

with  the  propagation  constant  defined  in  the  following  form  by  (35): 
k2  =  "K,  fC2(87),   where  K,=  lui^C    (88)f     ^  =  ^  (89) 

><  has  been  dropped  as  all  media     are  assumed  to  be  non-magnetic. 

Recalling  the  scalar  wave  equation  for  a  function  with  harmonic 

i  wt 
tine  dependence,    e  "     ,    we  have: 

(36)      V*U   +    K2U    =   O 

(31)      U--f(r)e 

The  spatial  part  f  then  satisfies: 

V2C(rW  K2-P(c)  =   O  (90) 

The  Ilertz-Debye  potentials  are  constructed  such  that  they  are  solutions 
to  this  wave  equation.     Let  It  (r,  Bf  0  }  =  f(r)   represent  either 
potential,   7T,  orlf: 

v2tt  +  k^tt  so  (9D 

This  equation  is  the  one  that  must  be  solved;  in  spherical  coordinates: 

7(-^r,l)+7^^i  15  (         Je>rWsaytKir=o     (92) 

It  can  be  solved  using  the  technique  of  separation  of  variables,  where 
iTis  assumed  to  be  a  product  function: 

lYCr.}  =  R(r)®(e)§CeD  (93) 

■>  -&r  Crir)  =  ®  $  £  r  *  (94) 


=>  ^(SINS!s)  =  Mis^e!f)  (95) 

=>  far  -.  R  s  *£.  (?« 

(92)  beoones:^  g£  +_M- |^  (sm  .  |f )  ♦$§1.  ^  a  -K*R©J 

2 
Since  'fP  ^  0,    then  we  multiply  by  r  /ft: 

This  equation  is  equal  to   a  constant  which  may  be  checked   by  talcing  the 
derivatives  of  both  sides  with  respect   to  r,     ©  ,  <p  .     This  constant  has 
the  form  of  n(n  +■  1 ) : 

X|2  rR+KV=   n(n  +  l) 


x    ri!rPUKxrz  Q Jl /sw  e  dfi51  J IUTI    (97) 


Thus  the  radial  equation  is: 


tne  radial  equation  is: 


(93) 


"sing  (98),  equation  (97)  becomes: 

By  the  same  differentiation  argu-ment,  each  side  of  this  equation  is 


2 
equal  to  a  constant,  m 


(99) 


Thus,    the  Q  -angle  eouation  is:  __ 

Using  (93)   and   (99),    (97)   takes  the  form  of: 
Thus,    the    <p  -angular  equation  is: 

*g£2+m*£0O»O  (ico) 

In  the  above  equations,  n  may  take  on  any  integer  value  and  m  may  be 
any  one  of  the  following  set:  -n,  ...,  0,  ...,  n.  Equations  (99)  and 


(100)  combine  to  form  the  Spherical  Harmonics  which. are  defined: 


in  * 

where  A  is  a  normalization  constant  to  ensure  orthonormality  of  the 

n 

function  set  of  Y. 

The  solution  to  (TOO)  is  a  straightforward  combination  of  sines  and 

1   2 
cosines  as  m  is  always  positive: 

$mW=  amcos(m^)  +  (^  sin(mtf)  (l02) 

Equation  (99)  is  the  familiar  Associative  Legendre  equation  having 
the  regular  solution: 

®Ce)  =  Pn     (ms)  (,0,;, 

(m) 

where  the  P.    (cos  &  )   are  the  associative  legendre  polynomials. 

Using  Rodrigues'  formula  Lll  : 

Pn   (C65©)=(2n!)  (l-x1)   J7^^  (104) 

where  x  =  cos  Q  and  the  negative  values  of  n  are  included,  -n^m  ^  n. 
'•/hen  m  —  0,  then  (99)  is  the  legendre  equation  whose  solutions  are  the 
legendre  oolynomials: 

PBft)«av)-,gr)..fx-')n  •  (,o5) 

Equation  (104)  may  be  rewritten  as: 

P^ld-  0-x*>*faP«M  006) 

The  complete  spherical  harmonic  solution  includes  the  normalization 
constant  and  rewrites  the  azimuthal  solution  as  a  complex  exponential. 
These  functions  are  completely  or tho normal  as  well.  Thus: 

Y>,  4)  ,  L-  0  V^  &$  PrfVs  e) «  l'mi/  (,ov) 

The   (-1)     term  is  a  phase  factor  defined  as  the  Condon-Shortly  phase,    and 

may,  in  some  texts,  be  dropped. 

(:Iote:  in  the  remainder  of  this  paper,  I  will  use  i  in: 
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which  is  used  in  the  literature.) 

The  Radial  equation,  (98),  has  solutions  defined  in  terms  of  the 
Bessel  functions: 


%(*')*  \f\*r    J^v,CK^  (103) 

XnO)  =  "Vf^7    Mn+'/zCKr)  (10?) 

=>      rRn(r)--Cn^(KOHnXr(^)  (]10) 

These  are  the  Ricatti- Bessel  functions,  where  J  j  and  II  ,i  are  the  half 

'       n+2     n+-g 

integral  order  Bessel  and  Neumann  functions.  By  using  the  infinite  series 

definitions  of  the  Bessel  functions,  one  can  write  the  Ricatti- Bessel 


d 


^00  =  z""  & 


n 


sm  z. 


z 


Cm) 

jr*M  =  or*""(^/(^  (m) 

where  z  B  x  +  iy  (113) 

The  function: 


3„<Kr>  =   ftM**  X„Otr)*JZKrft^Clir^         (m) 

vanishes  in  the  limit  as  kr  approaches  infinity,  and  './ill  be  useful  in 

(2) 
later  work.  The  function  H  /(kr)  is  the  half  integral  order  Hankel 

n+2 

function  of  the  second  kind. 

Now  it  can  be  seen  that  the  solution  to  the  scalar  wave  equation  is 
characterized  by  the  numbers  n,  m: 

and   (102),    (103),    and   (110): 


=>  rirr\r,.»-  itJM'dtXnMttrScMdkffi&j     (lis) 


Hie  general  solution  may  now  be  obtained  as  a  linear  superposition 
tiese  particular  solutions. 

rVTCO  =.    1  1        rTT<m\r,e,cO  (116) 


nio   ms-n 
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Equation   (11 6)    is   the  general   solution  of  the  scalar  wave  equation  in 
spherical  coordinates. 

The  fields  in  and  about  the  particle  can  be  described  in  terns  of 
their  individual  pairs  of  potential  functions; 

a.  the  incident  wave  -    "fTj  3  K^ 

b.  the  interior  wave  -  Tf,  )  Tr  ^ 

c.  the  scattered  wave  -  1T,  5  u a 

The  sphere,   whose  origin  coincides  with  that  of  the  spherical  coordinate 
system  and  is  of  radius  a,   is  isotropic  and  homogeneous,   optically 
described  by  a  complex  propagation  constant,  k   : 

U7)        =>         Kl    =  mi  K0  (H7) 

The  isotropic,   homogeneous  medium  is  a  dielectric  described  by  a  real 
constant,   kp. 

Kt-   ™Lk0  (118) 

The  relative  refractive  index  is  defined: 

m-.  KJk%  --  mjmi  (119) 

The  incident  plane  wave  propagates  along  the  positive  z-axis  and  is  con- 
sidered to   be  of  unit  amplitude  with  polarization  axis  parallel  to   the 
x-axis, 

|£'i|  '  \e'Lkl    ]   *  |  (120) 

In  order  to     facilitate  matching  the  boundary  conditions  for  the  Debye 
potentials  on  the  surface  r  =  a,    (120)   is  written  in  the  form  of  (116), 
realizing  that  the  function    Xn(^pr)   becomes  infinite  at  the  origin 
where  the  field   (120)   does  not: 

r<=i\%y  ^TT)  VKzOP„(%se>«'*  (-21) 


^-irX^^S^O^^W 


(122) 
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The  unperturbed  incident  wave  is  now  well  defined.      Because  the  incident 
wave  propagates  along  a  coordinate  axis   (z),  parameter  m  =  1.      As  a 
result,   the  Debye  potentials  all  have  the  associated  legendre  poly- 
nomials of  degree  1.     The  internal  particle  fields  must  be  of  the  same 
form,   without  Xn(k,r),    and  must  also  be  finite  at  the  origin.      Also,   the 
arbitrary  coefficients  must  exist  in  order  to  properly  map  the  conditions 
at  the  boundary; 

rTT'r-"  i*  1 '""  7?&fi  tft  *„  U,r)  pA«»  •)«*  4      (123) 

r  TT{  r  IT?  "4,   l"-'     -S    dn  *«0<,r)  PnW(c*  ^   SIN  *        (12A) 

The  scattered  wave  must  vanish  at  infinity  without  entailing  an  infinite 
energy  and  the  function  nn(kr)   contributes  just  this  property  while 
maintaining  the  proper  form  for  IT (r).     Arbitrary  coefficients  are  a- 
gain  necessary: 

rTT'    =  "Kj   nir  L  7^77)&n  jn(Kzr)Pn     ( Cos  e)  cos  0  (125) 

The  boundary  conditions  at  the  interface  ensure  the  continuity  of  the 

tangential  components  of  E  and  K.  In  soherical  coordinates,  this  trans- 
lates  into: 

Ee    (r»a)  =  Ee    {r^o.)  (127) 

E&   (rza.)   =•  EqS    (r.-a)  (128) 

H^(r-.a)  =  h£}  (r--a.)  (129) 

H$  Cr=a)  =  H^(r--a)  (130) 

where  (l)  refers  to  the  sphere  interior  and  (2)  the  exterior.  The  mix- 
ture of  terms  in  TTj  and  TTj,  at  first  glance  makes  it  difficult  to  apply 
the  boundary  conditions.   However,  by  forming  the  appropriate  linear 
combinations  of  (31)  -  (36),  the  boundary  conditions  can  be  recast  into 


A2 


0  =  ~hUe  +h  (s"<eM 


d  =  ^£ 


-  *- (S,«  6  E £ 


Af0(r.O   -    At0Cr  =  a) 
CWCr:a)=    CC°Cr-<0 


CO 


.£>■) 


Dwfr--a)»    Dw(r;^ 

Note  that:  IT^  =    TT  r 

The  quantities  A  and  C  remove  r  Tf2  and  3  and  D  remove  rll(. 


131) 
132) 
133) 
134) 


an  equivalent  form  for  the  Debye  potentials.   Define: 

a=  ■—  (sihqEg)  +  t?  £^ 

dgf 

Then  the  boundary  conditions  (127)  -  (130)  become: 

135) 
136) 
137) 
138) 
139) 
140) 
The  decoupled 

equations  then  result  in  boundary  conditions  on  the  Debye  potentials. 

Substituting  (31 )  -  (36)  into  (131)  -  (134)  yields: 

**™°k[4?(rvil**«»fc[&r*$+hr9&[4-rlrrf    cud 

'B  -  «  shC&rtt&  sws  &&Mll+hi%r&r(rrr3]    (142) 
rG  .  MBJifrqrtrJ*  SiHeg;.  [tf,  rfjffci  ^.[yC.rlQ    (143) 

For  equation  (135)  -  (138)  to  be  valid  for  all  9  ,  0  on  the  coordinate 
surface,  r  =  a,  requires  the  various  bracketed  arguments  in  (141)  - 


(144)  to  be  equivalent  between  regions  (1)  and  (2). 


ThU{ 


(145) 
(146) 
(147) 

(143) 
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Now,   using  (139)   and  (140): 


JL 


vM*$rM*tf2 


br   LL"\J    -    dr    L^"i  -  »'i  sj  (U9) 


<V^;=K,(0r(tr;t^)  (151) 

Yff  r<   =  >f5?r(ffi*0  052) 

Equations  (14-9)  -  (152)  are  the  required  boundary  conditions  on  the 
Debye  potentials.  Substitute  (121),  (123),  and  (125)  into  (151): 

iff, ' "" « Cn  y«  M  ?"<-as 9)('"^ = f I, l  ""^  ft  M-»3«M  ^ e)t"  ' 

Since  all  terns  in  the  series  expansion  are  linearly  independent  of  each 
other  for  each  value  of  n  and  for  all  & ,  ft;   then  the  corresponding  terns 
in  the  series  must  be  equal: 


K2  XV1* 

The  quantity — ^ — »  can  be  reduced  by  using  equations  (37),  (89),  and, 

recalling  that  the  frequency  of  an  electromagnetic  wave  is  unchanged  in 
traversing  the  interface: 

k]_y£      -  K?  K?  X  ■'      X?.      .w      , 


K*    K,("    "    -  V<T  X"  ?  K'"     "    X?  =    <w    = 

Applying  similar  arguments  to  the  remaining  boundary  condition  yields  a 

K  i      P"\ .  K0  y 

term  of  the  form:   =  — *—=■  =  rrii/ro,  ss  rn 

The  four  boundary  conditions  then  become: 

C„y„'(Kla.)»  m^W)-^,,'^)]         (153) 

CnVn(K,a)=  [>"(K^)-an^rKz^]  (155) 

dnVn(n,a)  =£"#<*««■)- */iJ/»0<*aX|  m2        C156) 

This  set  of  equations  uniquely  determines  the  coefficients  a  ,  b   c 


UU 


and  d  .  However,  only  those  corresponding  to  the  scattered  uave  are 
important  here.  Define  the  following  parameters: 

~  =    K,a--  Ztrj,   =    21T^  (157) 

d  =.  Kta.  =  "i^a.  =  ivrmji. '  x  mmi  (,53) 

where  Aft  is   the  vacuum  wavelength. 

The  quantity  <x  is   called  the  dimensionless   size  parameter.      Solve 
for  c     from   (155)   and  use  the  result  in  (153).     Also  solve  for  d     from 
(156)   and  use  the  value  in  {l5^+)»     We  obtain: 

an  =    j/,  r«*  >  jvw  -•»  y%  s/3)  y(oT)  U5^ 

b    -     2*JLQ  td^d }  -  V"  £gj  %Lt£-  (160) 

n  "    m  J*  r-o  Pi?  ^  - y„  r^j 37  f«*J 

All  the  field  vectors  are  now  uniquely  described  given  the  known 

parameters  m  and<*.  This  data  enables  the  coefficients  a  ,  b  ,  c  ,  and 

*  *!   p    n. 

d^  to  be  computed  and  their  potentials  determined.  Finally,  the  vectors 

are  obtained  from  (31 )  -  (36).  The  formal  solution  is  complete. 

Practical  light  scattering  observations  are  performed  in  the  far  field, 

or  wave  sone.  This  condition  is  such  that  k„r  >J  n  where  n  is  the 

Ricatti-Bessel  function's  order  parameter.   In  this  approximation,  the 

previous  relations  become  simpler.  The  asymptotic  expansion  of  the 

Hankel  function  becomes  flj  : 

3n(Klr)*i  Cf»Oe.i«tr  (]6i) 

The  fields  become  transverse  due  to   the  radial   component  decay  as 

2 
(  X/r)      compared  to  the   (X/r)   dependence  for  the  non-radial  components. 

Finally,  using  (31)  -   (36)   with  (125),    (126)   and  (1o1),    (162)   for  the 
scattered  wave  only  (   the  % (h„r)   tends  to   zero  for  largo  k„r  with  re- 
spect to  Jfl(k  r)   as  does  r^   )   yields: 
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-CK,r 


EW-wft  .-i-fc'    w?S,(9) 


O63; 


e.="&«'37   ir<"**»^  064) 


where  the  amplitude  functions  are 
°°-       2  n  +  / 


S^e)  =  ?,,    S    fatM*  tn  ITnCtoS©^  (166) 


and  the  angular  functions  are: 


%\h  e  '  (16?) 

(168) 
he  energy  flow  in  the  scattered  wave  can  be  calculated  by  using 


rn(t«  e)  =  ^r  PBc0(«se) 


Poynting' s  theorem: 


S  =    ^(e«H*  -E^He)  (169) 

where  ';:  represents  the  complex  conjugate  of  the  quantity.   The  intensity 
function  is  defined  as  follows: 

i  (9)  =  J  SO)  J  (i7o) 

The  intensity  of  the  scattered  radiation  for  unit  incident  intensity 
polarised  in  the  0  and©  direction  are: 

J©  =1^fvT  COS  0  lz  (171) 


(172) 


Note  that  the  intensity  decays  as  1/r  as  it  should  for  this  spherically 
spreading  wave.   The  scattering  plane  is  defined  as  that  plane  which 
contains  the  incident  direction  and  the  direction  of  the  scattered  wave 
(©,00.   Then  i.  is  the  wave  perpendicular  to  the  plane  and  i„  is  the 
wave  parallel  to  it.   The  phase  relation  between  E©and  Z^being  arbi- 
trary is  the  condition  that  defines  an  elliptically  polarized  wave.   In- 
tuitively, one  can  think  of  each  polarised  component  coming  from  that 
direction  as  being  inherent  in  the  incident  wave.   Therefore,  the  Igis 
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caused  by  an  incident  wave  of  intensity  cos  p?  polarized  parallel  to 
the   scattering  plane. 

The  phase  difference  that  exists   between  these  two  parts  of  the 


ave  is: 

~r  Re(S,)Ws,)  -  Re  (Sj)  Im  (S,) 

tan*>    -     Re(Sl)Re(S2)tImCSl)ImrSa) 


(173) 


Unpolarized  incident  light  gives  rise  to  a  scattered  beam  of  intensity: 


8ir 

with  degree  of  polarisation  given  by: 


P  = 


ll  +  LZ  I 


175) 


Fixing  the  scattering  plane  in  space  to   be   the  yz-plane,    all  ob- 
servations by  definition  take  place  in  the   scattering  plane.      The  in- 
cident electric  vector  has  its  direction  in  the  xy-plane  at  an  angle  X 
to   the  y-azis.      The  x-axis   is  designated  vertical  and  the  yz-plane 
horizontal.     This  is   shown  in  figure  S.      Thus,    the  vertical   (V)   and 
horizontal   (H)    components  of  the  intensity  in  the  scattered  wave  for 
unit  incidence   become: 


Tne  Stokes  narameters   are: 


S0  =  7^1  (i,  S/a/zX  -h  cz  cos  2  X) 
S,  =  ^TT  £/  sw2X  "  U  ^z  *>) 
si    =  g^X-*  y/Th.      six  X-   cos  x   cos  S 


A\ 


S//V  X     COS  X    SiM  f 


(176) 
(177) 

(17S) 

(179) 
(130) 


°2,       ztrW 

Uhen  discussing  cross   sections,    one  is   speaking  of  the  amount  of 
area  that  the  particle  effectively  presents   to   the  incident  beam  by  its 
ability  to   spread  out  the  incident  energy.      The   complex  index  of 
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Scattering  plane  geometry 


scattered  ray 


Figure  8 


IS 


refraction  includes  absorption  as  well  as  scattering.  Define  the  ex- 
tinction cross  section  as  the  sua  of  the  cross  sections  for  scattering 
and  absorption  : 

^ex-t      =  &%ccl   f  ^ahs  (182) 

By  integrating  the  real  part  of  the  time  averaged  Poynting  vector  for 

the  total  external  field: 

S  =(E'tE!)X  (Hl,  HS)  (183) 

the  total  outward  flow  of  energy  is  obtained.  For  unit  incident  energy, 
the  various  energy  losses  represented  by  the  integral  resulting  from 
above  are  the  desired  cross  sections: 

"«**&%t  C2nt0  £K«IMUM  da*) 

=-eMt=T^"i  O+O  £Re(a^kn)3  (185) 

The  efficiency  factor  is  defined  as  the  cross  section  per  unit 

2 

geometrical  area  (ffa  )  of  the  particle: 

Terms  of  the  formjL.  X  rr*t   become  ^t.  since  «<  =  2TT<1|\  .     Thus: 

Qsc^r^T    (Zn*d  i\o.„\\lbn\ll  (,87) 

Q.^z^r'%,  O+O  {R*  t***bJ]  (138) 

Note  that  ^a**   and  o-pv^  are  independent  of  the  state  of  polarisation  of 
the  incident  wave. 

Kerker  Q2J  et  al.,  gives  an  excellent  review  of  the  difficulties 
involved  in  the  numerical  computation  of  the  previously  listed  scattering 

functions.   The  functions  needed  to  be  evaluated  are  a  ,  b  ,  tf  *  X  ■ 

n7  n'  un'  Jn 

Some  discussion  of  the  equations  is  given  in  the  Appendix. 

Using  approximations  for  small  arguments,  the  fourth  power  law 
attributed  to  Rayleigh  scattering  can  be  directly  obtained.   The  first 
three  scattering  coefficients  are: 
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a,  --Vzi  (■&&)  ~J  da?) 

*>i  =17s  i-  (w'-l)^4'  (19D 

plus  terms  in  higher  powers  of°*.  For  the  Rayleigh  criterion8*^  <1 , 


only  a1  contributes.  From  the  Appendix: 

at,  (^5  g)  -.  i 

T,  fcos  &")  =  cos  © 

(i65),  (166):     s,(e)=  -|  a, TT,  TcoS  ©) 

st(e>  |  a,  r,  Cc.se) 

(170):  £/C)  -I  K|z 

For  unpolarized  incident  light,   for  example,    sunlight; 
(174):    Iu    =fc"    §l«-.|tCl+^SX«) 

(139):    Xu   = 


(192) 
(193) 


W- 


ml-l 


OZThV 


i^ 


<i57>.X„Sf&»(xirf)fcJ^|V«,*»:> 
==>Iu-l?^|^|l(1  +  C0S'e)  (194) 

This  equation  accurately  predicts  the  blue  sky  phenomenon  as  the  shorter 
wavelengths  scatter  much  more  than  the  other  visible  wavelengths  in 
daylight.  Equation  (175)  yields  the  observed  degree  of  oolarization: 

O1-  u^  to^B 


P-- 
p-- 


la.,!1  +  i^i*  cos*  e 
-  cos'-e 


(195) 


I  +C05  *  © 

A  plot  of  this  function  on  a  polar  graph  indicates  a  maximum  of  1 , 

complete  polarization,  at  a  scattering  angle  of  90  •  This  has  long  been 
observed  in  the  day  sky.  The  error  associated  with  the  truncation  after 
n  =  1  is  5/o  if  a/X  is  kept  -^  0.05.  This  limit  is  the  required  Rayleigh 
criterion. 
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The  cross  sections  and  efficiencies  represented  by  equations  (184.)- 
(188)  give  the  total  integrated  differential  cross  section  with  the 
integral  formed  over  the  entire  solid  angle  sphere  about  the  particle. 

The  differential  cross  section  is: 

ASL  CO-)         n  lr\ 
d-Tl      -  (s  /  (  l  (196) 

where:    7      ~  number  of  photons  scattered  into  a  solid  angle  per  unit 
1  3 

time  per  particle 
V  .   -  number  of  photons  crossing  a  unit  area  per  unit  time 
d-fl.  -  Siae   de   dtf         9   the  differential  solid  angle 
The  intensity  is  obtained  by  multiplying  the  number  of  photons   at  fre- 
quency V  by  Planck's  constant,   h,    as  hv>  is  the  energy  per  photon. 

h   =  b.UU  X   10  Toule/W  (197) 

As  written,    equation  (196)  describes  a  scattered  energy  per  incident 
intensity.     The  scattered  energy  becomes  a  scattered  intensity  when 
divided  by  the  distance  to  the  observation  point,    squared,   for  then  the 
solid  angle  becomes   the  subtended  area. 

dA-.rzd^  (198) 

Utilizing  the  above, 

AT  — 

dX     "    c(A  (199) 

v/here  dl  represents  the  scattered  intensity  over  dA  per  unit  incident 
intensity.  The  equations  describing  such  an  intensity,  (171 )  and  (172), 
may  be  cast  in  the  form: 

I'ifcaW'^^  (200) 


where        g  (0)    =     COS  Z  tf    or  Sin  ^  tf 
t  (r*,<*,  ©)  =  j,  on  z 
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I  is  the  dl  per  unit  intensity 


Thus  * 

(201) 


This  equation,   when  integrated  over  all  angles  Q9$  yields  the  cross 
section  given  in  equations   (184.)   and   (1S5).      Hereafter,    the  symbol  &  will 
mean  dO/dAas  defined  by  (201).      This  is  in  keeping  with  the  notation  in 
the  literature.  t  \       /    \ 

cO.-.e.sO  -^prU«.-.«)jW  (2o2) 

This  area  represents  the  scattered  energy  per  particle  per  unit  intensity. 

The  theory  may  now  be  extended  to  multiple  particle  systems  with  the 
addition  of  the  following  definitions.  Let  n  be  the  total  number  of 
particles  in  a  region  of  differentiable  volume,  cross-sectional  area  and 
extent  given  by  dV,  dA,  and  dx  respectively. 
Then : 

dV  =  dxdA  (203) 

and  n  =  NdV  *  NdxdA  (204.) 

where   N  -  particle  number  density 

Assume  the  scattering  to  be  incoherent  among  the  particles.  Then 
the  total tf  for  all  the  particles  is  just 

^r   =  n  er  (205) 

where  o  -  scattered  energy  per  unit  intensity  for  one  particle 

The  amount  of  intensity  scattered  into  dA  is  the  scattered  energy 

divided  by  dA. 

q"  f)   _   scattered  intensity  (206) 

~5~A  "   unit  incident  intensity 

Another  useful  parameter  is  the  scattering  coefficient: 

tf  scattering  intensity (207) 

incident  intensity-path  length 
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where  the  path  length  ax  is  that  traversed  by  the  incident  radiation. 
Combining  the  above  equations: 

=>  (5    s  CrN  (203) 

Until  now,  r  represented  the  radius  of  the  radiation  field  and  a  the 

particle  radius.  Again,  to  be  consistent  with  the  literature,  r  must  be 

assigned  to  the  particle  radius.  Using  the  new  notation,  if  there  exists 

a  population  of  scattering  centers  each  described  by  II (r),  then  the  a- 

mount  of  scattering  from  all  the  particles  of  radius  r  between  r  and 

r  4-  dr  is: 

C^{3  =  cr  (m^X.nQ^)    dhl  (r)  (20?) 

The  quantity  dIT(r)  nay  be  written  as: 

dtitr)   =  ^  dr  (210) 

where  -~  is  the  particle  number  density  oer  oarticie  radius.   Then: 

O  T*  * 

d(3=  <r(^)dr  (211) 


Let  the  lower  and  uooer  radii  limits  be  r,  and  rn  respectively:  we 

1      2  J ' 

haYe:   /a    r'1  ,  ,    v!  [iW\j 

In  the  case  of  cf  =0,-^>  then/3  =  /3g„£  ,  the  extinction  coefficient 
for  a  volume  of  scattering  particles  defined  by  dV  and  equation  (212) 
would  represent  the  total  energy  loss  by  an  incident  beam  passing 

through  a  volume  of  depth  dx.      Assuming  an  x   (path  length)   dependence 

d-(r) 
for  —y1 — L  ,    then  the  amount  of  intensity  lost  per  unit  incident  inten- 
sity is:  i. 

S8     =     '*<*>**  (213) 

Tex)    =  X.  e*P   {-J   /3(s)di}  (2U) 
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where  I0is  the  incident  intensity  to  the  volume  at  x-0.     This  last  e- 
quation  is  known  as  Bouguer's  law. 0>3J         It  describes  the  degradation  of 
the  incident  bean  as  it  traverses  a  material  path.      In  order  to  express 
the  scattered  intensity  into  a  solid  angle  di^we  recall  that  the  in- 
tensity scattered  per  particle  into  solid  angle  &£l  was  given  by  equation 
(200).     Let  the  parameter  r  in  that  equation  be  replaced  by  R  as  r  now 
describes  the  particle  radius. 

z^-feq?  IC m.*,«) 3 W  (2,5) 

As  before,  the  number  of  particles  in  the  cell,  II,  applied  to  equation 
(215)  yields  the  total  scattered  intensity  per  unit  radii  per  volume; 
then,  for  all  particles  of  radius  r  between  r.  and  r~: 

=>  I  =  1(1*1,0,00 

For  unpolarized  radiation: 

i  (no,  «*,  o)3(#)   s  £  (it  f  lz)  (217) 

where  i  and  i«  are  given  by  equations  (165)  -  (163)  and  (170).  For  the 
polarized  radiation: 

Um,<*  ,*)$&)  =  h(i-\toH%  # +it  cos  V^       (21S) 

Equations  (212)  and  (216)  are  the  specific  results  of  the  theory 
used  to  investigate  aerosol  parameters  and  will  be  discussed  in  the  next 
two  sections  of  this  paper.  First,  some  comments  on  the  collective 
scattering  phenomenon  are  in  order.  It  is  well  known  that  the  more 
ordered  the  scattering  centers  become,  the  more  coherent  the  scattering 
due  to  the  phase  relations  between  sites.  [Ts]   This  leads  to  a  reduced 
backscatter  with  corresponding  increase  in  forward  scattering.   In 
crystalline  substances,  there  is  essentially  no  backscatter,  when  off- 
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resonance;  and  the  incident  bean  merely  undergoes  propagation  velocity- 
changes.  The  condition  met  by  a  gas  at  atmospheric  pressure  or  less, 
where  the  number  density  is  low  on  a  radiation  wavelength  scale,  ensures 
that  the  scattering  is  incoherent,  allowing  the  summing  of  intensities 
from  all  the  scattering  centers.  The  Mie  theory  for  a  collection  of 
particles  requires  a  random  distribution  of  scattering  sites  to  achieve 
this  incoherent  intensity  accumulation.   Quantitatively,  the  near-field 
interactions  between  the  particle  may  be  ingored  if  the  average  distance 
between  sites  is  at  least  ten  times  their  radius.  This  criterion  is 
satisfied  even  for  very  dense  fogs  and  the  independent  scattering  model 
computations  are  valid  for  practical  tropo spheric  scattering.  However, 
significant  number  densities  will  lead  to  a  large  flux  of  previously 
scattered  light  in  the  region  and  add  multiple  scattering  effects  to  the 
received  signal,  A  quick  check  on  the  validity  of  the  single  scatterer 
model  is  to  verify  linearity  between  the  observed  effects,  I  and/3,  and 
the  number  density  of  scattering  particles.   If  the  relation  is  non- 
linear, one  must  suspect  the  existence  of  coherent  and/cr  multiple 
scattering.  Alternatively,  if  the  relationship  is  linear,  it  is  certain 
that  the  independent,  single  scattering  conditions  are  present.  A 
measure  of  this  condition  may  be  obtained  from  equation  (214.) .  (24]    Let: 

Tx  =  -J[    /3(s)ds  (219) 

be  the  optical  deoth  of  the  medium,  then: 

160  =  £0e  (22°) 

The  conditions  are: 

Tx  <  O.I   (transmission  >  9055)  -  single  scattering  predominates 
O.l -^  Tx  C  0.3   (transmission  75  to  95%)   -  multiple  scattering 

effects  may  be  removed  via  correction  to  the 
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single  scatter  model 
Qj^  T**   (transmission  <15%)  -   the  complete,  non-ana- 
lytic, approximate  multiple  scattering  theory 
applies 
Virtually  all  aerosols  in  the  troposphere  satisfy  the  first  criterion. 
The  literature  abounds  with  material  containing  graphs  of  the  Mie 
scattering  functions  for  many  specific  cases.  The  Appendix  outlines  one 
such  case  for  a  specified  set  of  conditions.   However,  some  general 
information  can  be  reduced  from  this  plethora  of  data.   The  amount  of 
forward  scattering  is  about  equal  to  that  of  the  back scattering  at  a 
size  parameter  («*")  equal  to  0.01.  As  ©*  is  increased,  the  amount  of  for- 
ward scattering  also  increases  with  a  like  decrease  in  the  backs catter. 
At«<-  1,  the  forward  scatter  is  on  the  order  of  100  times  that  of  the 
backscatter.  Increasing  »*  still  further,  there  exists  a  unique  scattering 
maximum  with  the  existence  of  smaller  subsequent  maxima  and  minima. 
The  scattering  function  increases  smoothly  as»*^  (Rayleigh)  up  until  o£* 
a/0.5»  Thereafter,  the  structure  is  highly  dependent  on  the  complex 
index  of  refraction.  The  large  oscillations  in  the  behavior  as»£  is 
varied  become  smoothed  as  the  absorption  index  increases.   As  the  num- 
ber density,  as  a  function  of  the  particle  radius,  takes  on  more  of  a 
detailed  structure,  there  is  a  further  smoothing  in  the  scattering 
functions  as  small  differences  tend  to  be  averaged  out.   Conceptually, 
one  has  at  hand  either  an  actual  particle  distribution  or  a  specified 
model  of  such. 

Given  rl(r),  for  each  r,  we  can  compute  the  required  I  (**,©)  and  sum 
the  total  intensities  using  IT(r)  as  a  weighting  distribution.  This 
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total  1(^,0)  then  appears  as  the  smoothed  function. 

Like  the  scattering  function  publications,  there  are  also  a  large 
quantity  of  analytical  number  distribution  models  in  use.   The  verifi- 
cation and  application  of  those  models  has  been  the  primary  purpose  of 
much  of  the  aerosol  experiments  conducted  in  the  past  20  years,  and  the 
work  is  continuing  today.   As  computers  become  faster,  real  time  inver- 
sion analysis  './ill  lead  to  a  more  accurate  determination  of  the  ritrue'r 
number  distributions.  Two  such  models  shall  be  mentioned  here.  Some 
natural  aerosols  may  be  described  with  a  two-parameter  function  utilizing 
a  multiplicative  scaling  coefficient,  c,  and  a  shaping  parameter ,1),  as 
presented  by  Junge     :  %  . 

KrJ--J-T^Cr  (22!) 

•..here   """,'  '  ^ '    is   as  previously  defined  in   (210)   ana  V  may  vary  from  about 
ci  r 

2.2  to  A.O.      As   an  example,    for  silicates  in  trie  O.l^m  to   10,0  uv.  range, 

one  may  use: 

r,   =    O.O^  yum  ,      rz  -    IO.O/*/n   ;     v   =3 

Another  such  model  is  in  reality  a  set  of  models  that  includes  the  Junge 
model  as  a  special  case.  Deirmendjian  UJ  has  proposed: 

tMmmfi  =  ftf -.-<*•*        '  (222) 

(note:   the  <=*"  used  here  is  different  from  the  size  parameter.)   The 

versatility  of  this  approach  is  the  presence  of  four  adjustable  parameters 

(a9oe ,   b,  Y)   allowing  for  a  wide  range  of  experimental  data  to  be  fitted. 

Three  specific  cases  of  this  model  are: 

Haze   G   (continental):        t  **  T 

-b\Tr~ 
Haze  M   (maritime)        :       r  ~  re 


Cloud 


r,       £    -1.5  r 
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'lost  of  the  other  proposed  model  distributions  utilize  some  Torn  of 
equation  (222)  as  their  basis  and  canned  routines  for  their  computer 
modelling  are  readily  available. 

Doppler  shifts  and  Brillouin  scattering.   The  presence  of  random 
(thermal)  acoustic  waves  in  an  atmospheric  cell  will  add  a  degree  of 
order  whereby  Brillouin  scattering  of  the  incident  wave  can  occur  off 
of  the  resultant  density  fluctuations.  However,  this  effect  is  extremely 
'.real:  for  gaseous  materials  and  the  degree  of  coherence  added  to  the 
structure  of  the  gas  is  generally  undetectable.  Similarly,  this  motion 
of  the  gas  will  cause  Doppler  broadening  of  the  wavelength  of  the  in- 
cident radiation  adding  to  the  natural  line  width  of  the  source.  Even 
for  the  sharpest  sources,  lasers,  the  natural  linewidths  are  on  the 
order  of  several  MHz,  Whereas,  for  almost  all  atmospheric  conditions, 
the  expected  Doppler  shifts  are  on  the  order  of  a  few  tenths  of  a  MHz, 
Therefore,  the  shift  is  generally  lost  in  the  natural  iinewidth  and  the 
relative  motion  between  the  detector  and  the  atmospheric  cell  may  be 
neglected. 

III.  LIDAR 

The  sophistication  of  the  laser  systems  available  today,  coupled 
with  the  fast  computer  processing  that  can  cheaply  be  obtained,  has 
singled  out  the  laser  radar,  LIDAR,  system  for  increased  use  as  a  tropo- 
spheric  experimental  device.  The  LIDAR  should  have  the  ability  to  per- 
form a  wavelength  scan  with  selectable  bandwidth  receivers.  It  is  the 
atmospheric  attenuation  of  the  laser  pulse  and  the  introduction  of 
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random  phase  and  amplitude  changes  that  combine  to   form  the  scattered 
pulse,   or  echo.     The  subsequent  inversion  of  this  data  then  allows  the 
experimenter  to   assess  the  optical   characteristics  of  the  atmosphere 
on  a  real-time  basis.     With  the  system  geared  to  perform  azimuthal  and 
elevation  scanning  as  veil  as  range  gating,   one   can  then  map  out  the 
atmospheric  optical  behaviour  within  the  visible  hemisphere.     The  two 
most  important  atmospheric  parameters  to  be  considered  for  a  pulsed 
LIDAR  system  are  the  extinction  and  scattering  coefficients.      The  total 
extinction  coefficient  may  be  written: 

A  =    aR+fl*+flAt/36  (223) 

where       R  -  Rayleigh  (molecular)   scattering 

M  -  Hie   (aerosol  +  other)    scattering 

A  -  Molecular  and  Aerosol  absorption 

0  -  Other  processes   (Raman,   fluorescence,    etc) 
The  other  processes  taken  together  are  typically  two  or  more  orders  of 
magnitude  weaker  than  the  Rayleigh  effect  and  tend  to  get  lost  in  the 
signal  noise. 

Rayleigh  scattering  is  obtained  from  the  Hie  solution  assuming  e>c<<l ; 


using  eGuations   (184)   and  (139)  -   (191): 
*.    r       \        V"         (  M  "\   \  m*--i   I  1 


^TTj  <*h  (224) 

Rayleigh  scattering  for  air  molecules  assumes  m  is  real  and  m  #1 ; 
*  =  £1T«/jl: 


8 


jp  (225) 


where  V  =  (4/3)fTa  ,  the  volume  of  one  molecule.  Let  N  -  1/V,  the 

s 

molecular  number  density,  then: 


3K 


^fi UJjl0,^  (226) 
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where  u  =  (6+3£)/(6-7^),  the  depolarization  factor,  required     (227) 
to  account  for  the  anisotropy  of  air  molecules,  for  which  cf  -  0.035.  [23] 
The  quantity 

g(©,0O   *  o.0b07i  Ci  +  C932cos2e)  (228) 

also  accounts  for  the  anisotropy.   Let  {°c^  be  the  air  molecule  number 
density  (vice  II),  then: 

/3r  =■  /^^,Cm,e)  (229) 

where  it  is  assumed  all  air  molecules  have  approximately  the  same  radius 
of  10"^i»n. 

The  Mie  extinction  coefficient  is  dependent  on  particle  radius 
through  equation  (212)  whose  only  unknown  is  dll(r)/dr.  The  unique  sharp 
wavelength  feature  of  the  LIDAR  allows  one  to  work  well  away  from  the 
characteristic  absorption  lines  of  the  sample.   Conversely,  it  also 
allows  one  to  type  the  scatterer  composition.  For  example: 

Laser  1:   \  on  an  absorption  (resonance)  line  of  a  particular 
species; 

=>  A3,  =  /3Rl  +  firm  *  @m  +  @f\  (230) 

where  p^  is  the  resonant  absorption  extinction  coefficient  which  is 
strongly  dependent  on  \.  with  tabulated  cV.  and 

The  (°   is  the  number  density  of  the  species. 
Laser  2:   a.. just  off  the  resonance; 

(3Z  --   &Rl  t  &mz   +  (3^  («2) 

Due  to  the  relatively  weak  ^  dependence  for  &$_^  @r<\\  WcL   P^ 

0,-Gt    a.  /Z'A  (233) 

For  an  analysis  cell  a  distance  R  away: 
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(214) 


i(i0.ioe-C'8Wds 


and  back  at  the  receiver:   « 

Then:   — r  -  q^~  (235) 

Thus,  the  (°'  for  the  particular  species  may,  in  principle,  be  found  from 
this  equation  using  the  inversion  techniques  outlined  in  the  next  section. 
For  the  typical  magnitudes  of  the  considered  cross  sections,  see  Table  II. 

TABLE  II 

Process 
Mie 

Rayleigh 
Raman 

Resonance  Raman 
Fluorescence 

Similar  to  radar  in  operation,  the  LIDAR  has  a  performance  equation 
whereby  the  system  characteristics  and  the  appropriate  atmospheric  para- 
meters wheih  contribute  to  the  signal  are  related. [23]   Consider  a 
pulsed  LIDAR  operating  in  monostaticmQ.de  (receiver/ transmitter  co-located) 
being  used  to  analyze  a  specific  atmospheric  cell  via  its  optical  prop- 
erties.  The  problem  geometry  is  outlined  in  figure  9.   In  this  figure: 

h.  -  cell  height 

R  -  slant  range  to   cell 

n   -receiver  solid  angle 
r 

A   -   transmitter  solid  angle 
ij 

L  -   cell  length  "gated"  via  receiver  on-time 
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hemostatic  laser  radar 
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Figure  9 
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A  -  cell  area  seen  by  receiver 

r 

a 


A,    -  cell  area  seen  by  transmitter 
The  following  derivation  shall  consist  of  counting  photons  at  each 
stage  along  the  beam  path;    converting  to   intensities  has  been  previously 
discussed. 

1,  fl,  -  number  of  photons   transmitted  oer  oulse 
•^L  -  photon  frequency  (transmitted) 

hVt.  -   energy  per  photon 
N.hV,    -   energy  per  pulse 
P.   -  laser  outout  power 
£*x  -  pulse  width 

=>     N,     a     S£  (236) 

2.  For  monostatic  LIDAR,    trie   scattering  angle  is   1S0  degrees    (ff  ) 
for  backsca.tter.      Assume  that  all  bachscattered  photons   are  received: 

=>  XI R    ^     nt  (237) 

The  ceil  volume  is  determined  by  the  LIDAR  timing  sequence  and  is: 

^=Lflt  (233) 

The  solid  angle:  n   *    f\/R2  (239) 

V  =  LiltR2  (240) 

The  pulse  width  constrains  the  maximum  system  resolution  by  limiting  the 
minimum  cell  length  achievable: 

Lm.n  =  ifeCcTj  (241) 

where  c  is  the  speed  of  light  in  the  medium. 

This  condition  assumes  that  the  pulse  tail  reflecting  off  the  front  face 
of  the  cell  does  not  interfere  with  the  pulse  front  reflecting  off  the 
rear  of  the  cell. 


3.  -•-  -  number  of  photons  incident  on  the  ceii 

(3a.  -  atmospheric  extinction  coefficient  for  the  given  \,  (")),) 

transmitted 

NlL  =Nte"^  ^tW^5  (242) 

However,  to  account  for  beam  forming  requirements,  assume  an  optical 
system  efficiency  of  v?,: 

=>    NL-^Nte'^rs)as  (243) 

4.  The  backs cattering  coefficient  may  be  defined  in  accordance 
with  equation  (212): 

where  (°cr(TT)  is  the  usual  notation.   This  equation  represents  the  back- 
scattered  intensity  per  unit  incident  intensity  per  unit  solid  angle  per 
unit  path  length.   The  signal  returned  from  the  ceil  is  due  to  this 
process  alone.  3ouguersT  lav  ( 21 4-)  states,  for  1° a" (flj  assumed  constant 

throughout  cell  length: 

-  p  cr  (J7)  L 

X-I0e  (245) 

where  1  is  that  portion  of  the  incident  beam  that  passes  through  L  un- 
affected.     The  scattered  portion  is   I     =  1  -  I   : 

Hs-.HL(\-c'P<f(n)i)     S   °  (246) 

where  N  -  the  number  of  photons  backs cattered  and  a  ratio  of  intensities 
s 

is  the  same  as  the  ratio  of  the  number  of  photons.   Typical  tropospheric 
aerosols  yield: 

(O   C(T7)L  <<  I  (247) 

(if  (°cr(n)L  =  0.02,  overall  error  1''.) 

The  restriction  this  places  on  1  is  an  upper  limit  of  several  hundred 
meters;  L .   Thus: 
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=>  Ms"  =  MLt°cr(/r)L  ^9~> 

The  prime  on  H  indicates  the  number  of  photons  per  solid  angle  that  have 
been  backs cattered.   Combining  equations  (24-3)  and  (249): 

|V/S'  =  K/t  Y[tf^<r(rr)L  e*p  {-£ Gt  U)  ds  ]  (25c) 

5.  ITow  the  echo  must  pass  through  the  same  propagation  path  on  its 
return  to  the  receiver  as  it  did  on  its  way  to  the  cell,   however,  the 
extinction  coefficient  may  be  different,  (3}   due  to  the  accumulation  of 
processes  iihe  Raman  scattering,  Srillouin  scattering,  Doppler  shift, 
etc.,  though  individually  negligible. 

11     -  number  of  photons  at  receiver  mirror 

Nm*Ms'are"-/'kr««iJ  (251) 

(239),  (25c):  rRr  -1  ,  1 

rr\  ■ a i252J 

6.  The  optical   system  will   collect  the   echo  photons  received  at  an 
efficiency  of  7      and  convert  them  to  an  electrics!   signal  via  a  detector 
with  a  quantum  efficiency  of  (     : 

Mf  =    Umlr  7cj,  (253) 

where  N  -  number  of  .oho tons  received  and  counted  oer  transmitted  eulse, 
r  c         ' 


Mr=     Ht7tggW2^eJ£&ffigi^.  fa®*0'®]**} 

R5 


due  to  backs catter  from  the  atmospheric  cell  of  interest. 
Thus:  r     /Ri 

(254) 

This  equation  may  be  recast  in  the  form  of  intensities  if  the  photon 
population  returned  is  large  enough  that  individual  photons  need  not  be 
counted.  Let  D  be  the  optical  system  aperature  (Dj.  -  transmitter;  D  - 
receiver),  then  the  intensity  becomes: 
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Reciting  equation   (254-):  <■     (*r*  t  \    o  •  >"!    I 

Lr  -   ■ '     ■       j^a  \<5o) 

Thus,  the  important  paraneters  for  the  LI3AR  analysis  of  an  atmospheric 

ceil  are  the  scatterer  number  densities  and  associated  cross  sections. 

For  a  given  system  design,  one  has  available  the  quantities  I, , ^ ,  , 

(r>  (q  j   A  ,  D,,  D  ,  and  can  select  the  values  for  "v^ ,  L,  and  R.  If  one 

assumes  a  slow  moving  target  cell,  then  the  pulse  round  trip  time,  t  ,  is: 
tr  a  2R/c  (257) 

=>    fl=  OA)*r  (258) 

Thus,  one  determines  the  range,  R,  by  turning  the  receiver  optics  on  a 

time,  t  ,  after  the  oulse  is  emitted.   In  order  to  avoid  the  oroblem  of 
'  r' 

range  ambiguity,  like  radar,  one  must  ensure  that  t^  be  less  than  the 

period  between  pulses: 

tr  <  l/PRF  (259) 

where  PRF  -  pulse  repetition  frequency 
The  unambiguous  range  limit  may  be  written; 

K     =  (c/z)/PRF  (260) 

where  R  -  unambiguous  range 

u 

The  cell  length,  L,  is  determined  by  the  time  that  the  receiver  is  turned 
off  again  a  time  At  after  t  .   For  the  entire  cell  to  be  scanned,  the 
pulse  must  traverse  L  twice  in  At: 

1=  (c/2)4t  (261) 

and  the  receiver  must  be  turned  off  at  time: 

t  =  tr  +  A-t  (262) 

One  can  then  select  L  and  R  by  programming  a  timing  sequence  for  At  and 
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t  .  Finally,  the  interpulse  period,  T,  which  is  a  system  deadtime  be- 
tween pulses  available  for  processing,  becomes: 

T--  (i/PRF)  -  rt  (263) 

Now  one  measures  I  and  Vr  and  applies  equation  (256)  to  determine  the 
desired  quantities  (°<y  6t)  and  (3(s J.  This  inversion  is  covered  in  detail 
in  the  next  section.  Note  that  when  the  equation  for  tfO™^;  r)&>rf)   from 
(202)  is  used  to  replace  the  tftm,  A/,©,  $)  in  (244.),  then  equation  (256) 
for  I  has  the  correct  form  as  given  by  (21 6)  for  the  scattered 
intensities. 

Thus,  a  model  has  been  constructed  whereby  the  system  parameters  are 
related  to  the  cell  unknowns  in  a  definite  manner.  But  first,  some  con- 
sideration must  be  given  to  system  noise.  During  the  time  the  receiver 
gate  is  open,£t,  the  system  will  collect  and  count  the  background  sky 
photons  as  well  as  the  signal.  The  continuous  "dark  current"  in  the 
detector  yields  an  electrical  noise  as  a  result  of  dark  counts.  These 
counts  arise  from  photoelectrons  being  emitted  from  the  detector  photo- 
cathode  surface  even  without  an  illumination  present.   Each  of  these 
noise  sources  may  be  reduced  by  proper  design  of  the  equipment  and  ex- 
periment. In  the  case  of  a  detector,  its  dark  current,  I_,  is  either 
specified  or  may  be  measured  in  the  laboratory: 

XD  =■   NDeG  (264) 

where   iL  -  number  of  dark  counts  per  sec 

e  -  electron  charge  per  count 

G  -  detector  count  gain 
The  number  of  dark  counts  detected  as  part  of  the  return  signal  in  a 
time  At  is  NpAb,  and  : 
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ND  =  ID/(e&)  (265) 

Let  N„  represent  the  background  count  rate  such  that  N^At  is  the  num- 

D  D 

ber  of  background  counts  received  per  laser  pulse.  Assume  a  background 
spectral  brightness,  './^,  that  accounts  for  all  the  natural  and  artificial 
background  sources.  This  \I^  represents  a  power  density  per  solid  angle 
in  the  background  sky.  Then: 

fjB  =  (yv)  ikfkb  ax  (266) 

h  v 
where   v?  ^f)-nrjDr  -  previously  defined 

V  -  frequency  of  background  photons  (atX) 

A  A  -  receiver  bandwith,  centered  about  'X 

Let  II  be  the  total  number  of  noise  counts  per  laser  pulse;  using 
n 

equation  (261 )  for  At: 

Mn  =  (2L/c)(fVNB)  (267) 

Due  to  the  pulse  repitition  frequency,  PR?,  the  noise  count  rate: 

M=   (PRF)Nn  (26s) 

where   IT  -  noise  count  rate 

The  statistical  nature  of  noise  yields  an  effective  RMS  noise  count  rate 

of  VTT.  The  signal  has  a  count  rate  of: 

Sr  =  (PRF)Nr  (269) 

where       S     -  count  rate  from  signal  echo   alone 
r  to 

Thus,   the  total   counting  rate  is: 

S  =    pRF(Nr<-N)  (270) 

Talcing  the  RMS  of  the  signal  plus  noise,   one  can  form  the  signal-to- 
noise  ratio: 

S/N    =   5r/V3~  (271) 

Assuming  a.  one  second  time  averaging: 
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S/A/    *\fPRF    ^r/^TJI^  (272) 

If  the  noise  average,   N  ,    is  not  accurately  measured,    the  U     is  replaced 

by  211     (measured)   in  the  above  equation.     Note  that: 

1.  if   ^n=  0   =>    (s/A0o   =VCP«F;Nr  (273) 

2.  if    Nn=  Nr  ->  (S/N)L  =    Cs/A/)e//T  (274) 
Therefore,    the  S/1I  ~  ?RF  and  by  increasing  PRF,   one  increases  the  signal- 
to-noise  ratio.      However,   via  equation   (260),    the  unambiguous  rsnge  de- 
creases faster  than  the  increase  in  S/N.     Since  a  large  S/lI  is  desired 
with  good  range  resolution,   one   selects  a  detector  with  a  low  3L(lO   and 
designs  the  system  for  low  background  reception   (N-).      If  the  time  aver- 
aging    had  been  performed  over  a  period  of  T seconds   such  that  PR?  T 
pulses  had  been  transmitted  and  received: 

S/N  =  Vprft    Mr/v/7^T^T  (275) 

Combining  equations  (265)- (267)  -ana  (^7>;: 
_S       *    nTprft  A/r 

Examination  of  the  above  equation  and  (254-)  leads  one  to  some  genera] 
system  noise  qualifications. 

1.  One  desires  a  low  detector  I_  and  a  high^o. 

2.  Select  a  narrow  receiver  spectral  bandwidth, A  X;    another  reason 
for  the  choice  of  a  laser. 

3.  Reducing  L  reduces  the  noise  as  well,  but  the  transmitter  must 
put  out  greater  power  to  maintain  the  same  signal.  Then,  for  the  seme 


-j 


aver  at?  i 


e  power,  the  FRF  must  be  lowered  and  the  effect  in  S/H  is  only 
VPR?  as  compared  to  the  1  decrease  and  rTj.  increase,  increasing  S/IT. 
In  general,  high  power/low  PR?  yields  a  better  S/lT  than  a  low  power/ 
high  PR?. 
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4.  A  snail Slf  requires  a  lover  beam  divergence  at  the  transmitter 
and  an  optical  matching  between  the  receiver  and  the  echo  beam  pattern 
(another  laser  capability). 

5.  Without  any  noise,  there  still  exists  a  finite  (S/l-T)n  which 
represents  the  "shot  noise"  statistical  nature  of  the  signal.   It  must  be 
emphasised  that  oil  the  quantities  in  equation  (276)  represent  the  mean 
of  the  indicated  variables. 

6.  If      Mr   «^n      =>      S/N    ~Nr  (277) 
If       Mr    ^     J\|n      =>      S/M     ~  >fi£  (276 


jolecular)    and  hie   (aerosol)    scattering.      As   an  example; 

a.  single  water  droplet,    r  =   l^m   (hie) 

tfs   ~    !0"8  cm2 
(°    *•*     10     cm"3 
->   (3m  ~  lO"7^"1 

b.  nitrogen  and/or  oxygen  molecule,    (hayleigh) 

^s   ~    10        cm2 


7C-\ 


P    ~   10       am"3 

=>    @R  ~    I0"qcm' 


f|3o   to   the  total  (3 r    is 


on  "one  order  01   a  xew  oe] 


cent  and  any  quantitative  measurement  must  take  into   account  both  the 


■hiiT>T"n  A  -i  f. 


rat  .10 


of  the  total  scattering  Intensity  to  the  Rayleigh  scattering  intensity: 

TA  --  i-t/ir  (275) 

Some  researchers  attempt  the  measurement  of  this  quantity  independent  of 
the  intensity,  then  use  the  result  to  deduce  the  aerosol  I...   Another 
method  assumes  that  the  aerosol  concentration  above  a  certain  height 
(e.g.  ^-  30  Inn)  is  sufficiently  low  to  alio'.;  it  to  be  neglected.   Then 
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equation   (256): 

y.cxt,   one  computes   Ip  from   (256)   assuming  the  conditions  of  the  U.  3. 
Standard  Atmosphere  and  equation  (226).     This  is  the  expected  echo 
intensity  due  to  molecular  scattering  alone.      The  ratio  between  the  mea- 
sured  I_  and  the  expected  Ip  yields  a  scale  factor,    b,    that  corrects  the 
standard  atmosphere  for  the  local   conditions.      Using  this  factor,    one 
assumes   a  density  distribution  for  the  local   atmosphere  the  same  as  that 
for  the  standard  and  comoutes  the  exnected  value  of  I_   at  the   cell's 
location.      This  intensity  is  then  scaled  using  b.      The  measured  value  of 
IT  is  then  corrected  for  this  Ip(b)   and  the  result  becomes  L,,    the  Mie 
intensity: 

^M   =    ^-T   "    bIR  (2S1) 

For  other  means  of  separating  the  molecular  and  aerosol  intensities  from 
a  measurement,    the  reader  is  referred  to  the  work  by  Hall,    et.    al.   Qf] 
This  Rayleigh- atmosphere  assumption  is  the  starting  point  for  the  in- 
version process. 

The  Standard  Atmosphere  may  be  represented  via  the  following  equations; 
temperature  T(°F)  =  59.31  -  0.00357  (°f/k)  U  (H)  (2S2) 

pressure     P(a4m}=   1.022.  etp  (~H.  157  *  \o'5(  i/pOh  (ft))  (283) 

Mass  density     .  0  002M  15  e^p  ("3-3$  X  lCT50/ft)Kft±))  (284.) 

(slugs/ft^5)    -      ' 
where  h  may  vary  from  0  to  37,000  feet.{j3j        The  more  precise  form  for 

these  relations  are  also  represented  as  canned  computer  routines. 

The  critical  term  in  equation   (256)   is  the  P    Cf^fr)  where  the  sub- 
script M  refers  to   the  aerosol  distribution  alone.      The  extinction  com- 
ponent in  the  exponent  of  e  in  this  equation  represents  the  scattering 


71 


and  absorption  cross  sections  integrated  over  the  entire  4-Tf  steradian 
sphere.  If  one  assumes  some  representative  model  for  the  aerosol  dis- 
tribution, as  was  done  for  the  Rayleigh  case,  then  the  only  unknown  in 
(256)  is  i^j  t4L0y*  '-^e  ^w0  averaging  processes,  integrating  aver  4-TT  and 
integrating  along  the  path  R,  reduces  the  sensitivity  of  the  exponential 
term  in  (256)  on  the  actual  aerosol  distribution.  After  completing  the 
inversion  process,  one  may  input  the  best  estimate  for  this  distribution 
into  the  exponential  term  in  (256)  and  then  repeat  the  inversion.  How- 
ever, the  processing  time  involved  may  prove  to  be  the  final  constraint. 

The  LIDAR  experiment  yields  the  scattered  Mie  intensity  from  the 
aerosol  population  in  the  analyzed  atmospheric  cell  as  a  function  of 
the  pre-selected  parameters  at  the  laser  (e.g.  wavelength  and  polarization 
scanning) .  With  this  data  as  input,  the  system  is  prepared  to  perform 
the  inversion  routine.  Some  computational  needs  are  reduced  for  the  case 
of  laser  backscatter  as  the  scattering  angle  remains  fixed  at  180  degrees. 
This  backs cattering  condition  simplifies  several  calculations  as  out- 
lined in  the  Appendix. 

IV.  MATHEMATICAL  INVERSION 

The  assumption  made  throughout  this  paper  implies  that  the  experi- 
menter need  only  apply  the  exact  Mie  theory  to  his  observations  and  per- 
form a  simple  data  inversion  to  discern  the  aerosol  properties  of 
interest.  However,  this  inversion  is  by  no  means  simple.   Aerosol  prop- 
erties may  be  separated  into  two  distinct  catagories.  The  first  con- 
siders those  parameters  that  are  explicit  in  the  previously  defined 
scattering  functions;  i.e.,  spherical  aerosols,  aerosol  sizes,  low 
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number  densities,  and  the  relative  indices  of  refraction.  Then  there 
are  those  aerosol  characteristics  that  require  some  approximations  (in 
many  cases,  severe  ones)  to  the  theory;  i.e.,  non-spherical  shapes, 
anisotropy  of  the  shapes  with  respect  to  the  scattering  plane,  and  high 
number  densities  leading  to  multiple  scattering  effects.  The  first  few 
LIDAR  experiments  investigated  the  actual  aerosol  ranges  for  these  para- 
meters. They  found  that  in  excess  of  90%   of  the  real  aerosols  studied, 
the  spherically  symmetric  single  scattering  theory  was  valid  under  favor- 
able atmospheric  conditions.  These  conditions  require  the  absence  of 
dense  water  clouds  in  the  atmospheric  cell  and  a  depolarization  of  the 
incident  beam  limited  to  a  maximum  of  2%,   with  ],5%   of  less  actually 
measured. £23}   In  the  following  analysis,  these  favorable  conditions 
will  be  assumed,  as  will  the  requirements  of  single  particle  scattering. 
These  conditions  are  enhanced  by   the  common  field  of  view  shared  by  the 
LIDAR  transmitter  and  receiver. 

Recall  equation  (21 6)  for  the  scattered  intensity  a  distance  R  away 
from  a  cell  containing  an  aerosol  population  described  by  dil/dr  with 
radii  limits  r.  to  r,:    y  fAUt  \\ 

(216)    I  .t$5J  f  lc^,-<>©)S^(-dr7  dr 

The  first  simplification  is  to  replace  the  particle  radii  continuum  with 
a  set  of  distinct  particle  classes  each  with  a  characteristic  radius, 
number  density  and  complex  index  of  refraction.  Thus: 

th 
where   m.  -  j  *  particle  refractive  index 

r.  -  i   oarticle  radius 

3 

N.  -  j   particle  number  density 

J 
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A.  -  incident  radiation  wavelength 
©,  P   -  scattering  angles 
The  large  number  of  unknowns  that  need  to  be  determined  will  ensure  a 
non-unique  inverted  solution  from  the  limited  sample  of  observa.tions_.a- 
vailable.   There  is  also  the  added  difficulty  of  variance  in  the  sup- 
posedly known  data  whose  uncertainties  at  best  generate  unbounded  noise 
in  the  final  results.  Therefore,  some  attempt  must  be  made  at  reducing 
the  initially  large  number  of  unknowns. 

The  dependence  of  the  intensity  on  the  size  parameter,  has  been 
discussed  previously.  3riefly,  the  non-mono tonic  fluctuations  as  °<  is 
varied  is  the  distinctive  feature  of  the  Mie  theory  that  enables  the 
process  to  be  of  good  use  in  tropospheric  analysis.  These  maxima  and 
minima  are  reduced  by  increasing  the  field  of  view  and/or  the  transmitter 
bandwidth.   It  is  to  LIDAR' s  advantage,  then,  to  retain  this  complex 
structure  as  a  data  filter.  See  the  Appendix  for  a  sample  plot  of  Q,(oC\ 

The  complex  index  of  refraction,  m,  has  a  dispersion  relation  like 
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The  less  than  \%   change  in  the  index  of  refraction  is  typical  for 
dielectric  media.  The  change  ir  the  absorption, y ,  is  even  less.  Studif 
have  indicated  that  the  dispersion  in  m  is  balanced  by  the  randomness  of 
the  aerosol  sise  distribution  in  such  a  manner  as  to  almost  negate  this 
dispersion.  \23\ 

The  angular  dependence  is  a  pre-specified  set  of  conditions,  espec- 
ially for  the  case  of  a  LIDAR  backscatter  system.  This  system  consists 
of  a  linearly  polariaed  wave  undergoing  only  backscatter  observations. 

As  a  result  of  these  and  similar  considerations,  one  can  reduce  the 
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number  of  unknowns  drastically  by  making  only  a  few  simplifying 

assumptions.   These  have  also  been  verified  by  the  laboratory  analyses 
of  many  aerosol  samples,  natural  and  artificial. 

TA3LS  III 
Liquid  water  (@20  C)  index  dispersion;  m  =  n(l-iV) 

,0' 


in. 


4046.6  1.342724 

4471.5  1.339423 
5015.7  1.336363 
5460.7  1.334466 

5092.6  1.332983 
6562.3  1.331151 
7065.2  1.330019 

1.  The  measured  values  for  n  range  from   1.33  to    1.59   with  respect 

to   the  surrounding  air.      A  large   class  of  aerosols    (particulates)   may 
be   characterised  by: 

no  =  I.  5 4.0-  tO.OoO  (2:6; 

Thus,  assume  that  the  complex  index  of  refraction  is  both  known  and 
unique.  Only  the  particle's  radius  ana  number  density  is  allowed  to  vary. 

2.  Typical  non-precipitating  aerosols  in  the  troposphere  are  described 
with  their  radii  between  0.01 ^n  and  15^  m.  Their  number  densities  can 

be  represented  by  distributions  containing  only  a  few  parameters.   The 
important  'model  distributions  were  presented  eajrlier  and  the  usual  goal 
of  an  inversion  is  to  best  fit  one  of  them.  The  initial  estimates, 
where  required,  of  particulate  properties  usually  includes  such  models. 
The  set  of  unknowns  reduces  to  the  earticle  number  densities  and 
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their  radii.   By  the  use  of  pre-selected  values  for  r,  one  only  has  to 
solve  for  the  distribution  functions  granting  that  the  radii  range  is 
properly  considered.  Therefore  the  parameter  of  primary  interest  is  the 
scatterer  number  density. 

Define  the  Mie  function,  K: 

Kfm,r13LJe,rf)s^pjj  i(mt«,e)  g(et)  (287) 

and  the  particle  size  distribution  function,  f(r): 

}  =   ~7r —  (2S3) 

The  scattering  intensity,  equation  (21 6),  becomes: 

100  =  J  \Ur)  fCO  cJr  (289) 

where  x  represents  the  set  of  user  defined  parameters  (n,  \,Q,  0).  The 

value  l(x)  is  measured,  the  quantity  K(x,r)  is  computed  and  the  unknown 
f  (r)  must  be  solved  for.  This  relation  is  classified  as  a  Iredholm  in- 
tegral equation  of  the  first  kind.   The  kernel  of  the  equation  is  K(x,r) 
and  l(x)  represents  the  transform  of  f(r)  utilizing  the  stated  kernel. 

There  exists  only  a  finite  set  of  measurements,  l(x),  over  which  the 

th 

inversion  may  be  performed.  Let  the  subscript  i  represent  the  i 

measurement  of  a  set  of  n  observations: 

I(*0  =  9(*0«9t  =  J       K(^,r)f(r)  dr  (290) 

The  notation,  g.,  has  been  introduced  and  will  be  used  later  in  the  paper* 
The  inversion  process  is  by  nature  an  unstable  one  and  this  may  easily  be 
demonstrated  •  [22J    Define: 

SeOuO=    fKttoOpM+CswFr]    dr  (291) 

Jo 

where       i  =  1,    2,    ...,   n  F  -   1,    2,   3,    ... 

and  C  is  an  arbitrary  constant.   Letting  r  =  0  and  r?  =  IT  in  equation 
(290): 


76 


5F(*i>f)    =3^x^+    C  JoSLn(Fr)k(tL,r)dr  (292) 

Since  the  sin(Fr)   acts  as  a  weighting  function  for  the  kernel  on  the 
interval   (0,11),    as  7  increases  without  bound,    the  second  integral   tends 
to  0.     Thus: 

SF(*ti*)P-..     "9  60  (293) 

One  can  approximate  the  transfomed  values  of  f  with  the  function  S^  as 

close  as  possible  by  talcing  F  large  enough.      Yet  the    jOsihTr     function 
may  be  made  arbitrarily  large  and 

Hr)    =  f(r)   +  CslnFr  (294) 

will  transform  as  f(r)  with  little  to  no  resemblance  to  f (r)  for  large 

values  of  G.  Therefore,  one  can  generate  transforms  that  come  as  close 
to  the  desired  values,  g(x),  within  any  desired  accuracy  as  one  may  wish 
and  still  have  no  idea  what  the  accuracy  with  respect  to  f(r)  may  be. 

The  technique  used  to  remove  this  instability  is  to  require  the  ad- 
dition of  a  constraining  relation  on  f(r)  other  than  equation  (290). 
This  may  be  done  either  explicitly  or  implicitly.  The  proper  relation 
must  contain  the  desired  qualities  of  f(r)  requiring  a  high  degree  of 
a  priori  knowledge.  Out  of  the  near- infinite  set  of  possible  solutions 
for  f (r)  generated  by  the  inversion,  the  constraining  relation  will  act 
as  a  filter  to  obtain  the  unique,  or  near-unique,  solution. 

The  presence  of  measurement  error  in  the  g(x.)>  larger  than  some 
quadrature  error,  allows  the  use  of  an  appropriate  quadrature  approx- 
imation of  the  form: 

for  k  intervals  of  the  argument  y  and  the  w.  are  the  quadrature 
weighting  coefficients  with  increments  Ay,- •  For  equation  (290): 
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3 (*<>=  %  Wi,Kf*i)'3)f(rOArj  (296) 

Oil  J 

where  i  represents  a  measurement  and  j    a  radii  interval. 
Define;    for  n  measurements: 

1.  column  vector  G:     g.   =  g(x. )       i.  -   1,    . ..,   n;  (297) 

2.  coefficient  matrix  A:     a. .  a  w.  .K(x.«r.)  (29S) 
j  —   1,...,  1c;     k       radii  intervals 

3.  column  vector  N:     N.  =  f(r.)  Ar.  (299) 
Then,   equation   (296)  may  be  written  as  a  matrix  multiplication: 

G  =  AIT  (300) 

The  value  N.   represents  the  total  number  of  particles  per  unit  volume 
with  radii  between  r.    and  r.  +   Ar.,    as  can  be  seen  via  (233): 

Nj   *J       f(r)3r  =  f(rpRrJ+Arj)-ri]  (301) 

and  f(r. )   is   considered  constant  over  the  interval  j.      It  is  further 
assumed  that  the  scattering  functions  have  only  small  variations  with- 
in each  sub- interval  radii  range. 

The  problem  has  now  been  reduced  to  solving  for  the  k  unknowns  in 
N  given  the   (k)x(n)    computed  values  in  A  and  the  n  measurements  in  G. 
It  must  be  emphasised  that  a  quadrature  approximation  will  introduce 
serious  error  in  f(r)   unless  the  quadrature  weighting  formula  is  of  an 
accuracy  higher  than  that  of  the  expected  measurement  error.     Let  nn 

be  the  number  of  independent  observations  made  such  that  n^  <   n.      If 

u 

k  -  nn,   then  II  may  be  solved  for  exactly: 

(require  A  to  be  non-singular) 

=>    a"1  A   =    A  A"1  =  I  (302) 

whercx   1  in  a  matrix  equation  represents  the  identity  matrix, 

(300):   f\H  -   G   =>  NJ  =  A"'G  ,   (K^nJ  (303) 
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A 

where  M  represents  the  interpolated  solution  vector.   However,  if  k 
K  n~,  then  II  may  be  solved  by  minimizing  the  square  of  the  error  (least- 
squares)  : 

(require  A  A  be  non-singular  =y   A  :  transpose  of  A) 
(300):    (flrA)N  =  AT£r  (304) 

this  step  forms  a  square  matrix  which  may  be  inverted, 

M  =  (ATA)"'Ar&  ,     (KirO  (305) 

That  this  last  equation  actually  minimizes  the  error  may  be  visualized 
by  considering  an  appropriate  spatial  geometry  which  contains  the  vector 
G  and  the  vector  AIL  where  IL  ^  N.  See  figure  10.   The  vector  G-AIL  i 


s 


seen  to  be  representative  of  the  error  vector  between  G  and  All  .   The 
minimum  error  occurs  when  this  error  vector  has  the  smallest  magnitude, 
thus : 

CAN,)  ■(&-/!«.)  =0  (306) 

(the  •  represents  the  vector  dot  product). 

A 

Then  for  any  All  in  the  range  of  A,    the  vector  AIL    is   still  perpendicular 
to   it: 

W)-(G-AN)  =   O  (307) 

Writing  the  dot  product  as  a  matrix  multiplication  requires: 

(AN,)T(G-flW)=  o  (308) 

=>  h),T(flrG-ArflN)  =  O  (309) 

And  NT  ^  0  since  it  represents  any  TL   vector  in  the  range,    thus: 

(304)  f\Tf\N    =■  f\7Qr 

In  the  work  that  follows,    when  writing  the  inverse  matrix  A     ,    if  the 
matrix  is  not  square,    the  quantity  (ATA)      A     is  implied. 

The  procedure  as  indicated  by  equations    (303)    and   (305)   are  in  large 
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n-Space  vector  geometry 


range  01 
operator  A 


Figure   10 


unsuccessful  as  they  do  not  contain  the  constraining  criterion.  Some 
further  problems  concern  the  size  of  k. 

A 

1.  Small  k:   the  vector  II  does  not  accurately  represent  N  due 
to  the  r.  set  sparsity  and/or  the  crudeness  of  the  quadrature. 

A 

2.  Large  k:   the  components  of  N  take  on  an  oscillatory  be- 
haviour irrespective  of  N  and  the  coefficients  in  A  become  particularly 
sensitive  to  round  off  error. 

This  sensitivity  in  coefficients  is  -est  characterized  by  the  con- 
dition number  of  A  (or  any  matrix)  £5}  : 

concf(A)  =  |fl|  JA~'|  (3io) 

'/here        |A|    -  norm  of  a  matrix 

The  norm  of  a  matrix  may  be  defined  in  any  one  of  three  ways  depending 

upon  the  desire  of  the  user  to  highlight  specific  aspects: 

I  o  I     -    nncu         %     ,         I  ,_„,, 

'        ISJ^K    tVi    I  CL c j I     :    maximum    ctiumtj  -  sum  1311; 

v/here  a.  .)    -  magnitude  of  the   component  a.  . 

,  „i      rr\ojC         4  1  ^   I  .  .      „. 

in  loo  =  l^i^n  j>,  l^jl  :  maximum  ro-.^sum  (312) 

(normally  used)  |fl|e  =■  (  v  S  I  &i^\   J  :  Euclidean  norm  {313) 

Because  A   must  be  computed  and  will  contain  at  least  the  inaccuracies 
in  A,  the  cond(A)  will  be  costly  to  compute  and  inaccurate  as  well. 
However,  the  test  is  the  magnitude  of  the  cond(A).   If  large,  the  matrix 
A  is  termed  poorly- conditioned  and  the  likelihood  of  a  successful  inver- 
sion decreases  rapidly  to  zero.   Corrections,  like  scaling  or  partial 
pivoting,  must  then  be  applied  to  A.  Hote,  for  the  identity  matrix, 
cond(l)  =  1; 

cond  (A)  >y   I  (3U) 

since  A~  must  exist.  Other  properties  of  the  norm  arc: 


IAl£0  ^^ 

\(\\  =    O     'iff     Pl*Q  (316) 

|/3fl  |    =  A  \f\\  i    &  >  O      dnd  a.   number-  (317) 

IA+BI  ^   Ifll  +lBl  (3i3) 

iBfll  <  I6|  IAI  (319) 

The  ccnd(A)  is  most  useful  in  bounding  the  error  expected  in  the  computed 
solution  and  relating  this  error  to  the  magnitude  of  the  residuals. 
Define: 


AX  =  R 


U2CJ 


where   Z  -  exact  solution  vector 

3  -  approximate  solution  vector 

3  -  vector  of  known  quantities 

A  -  operator  matrix 
The  residuals  are: 

R  =  3  -  AX  (321) 

'■here   R  -  residual  vector 

A 

culu .  _j  —  .I.  —  a.  \j(LtL) 

Then:  R  =  A3  (323) 

using  3  by  inverting  equation  (323)  and  the  norm  relations  (315)-(319), 
one  may  incorporate  the  cond(A)  as  an  error  measure: 

A 

Therefore,  the  relative  error  in  X,  contained  in  3,  can  be  as  great  as 
the  relative  residual  times  the  cond(A)  and  as  small  as  this  same  re- 
sidual divided  by  cond(A).   Tor  large  cond(A),  the  residual  yields  iibt 

A 

to  no  information  on  the  accuracy  in  X.  Equations  (320)- (323)  may  bo 
used  to  define  an  iterative  improvement  algorithm  (I1A). 
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A  A  i 

compute  X  oy  X  -  A     3  1 325) 

A 

(321 )        2.      compute  R  by  R  =  3  -  AX 

3.      compute  E  by  E  =  A_1R  (326) 

4..      correct  X  by  X  =  3  +  X  (327) 

If  the  system  matrix  A  has   errors   in  its   coefficients,    the  system 
may  then  become  ill-conditioned  even  if  the  residuals   are  small. 

A  =  A  +  H  (328) 

".;here       A  -  true  coefficients 

H  -  error  in  coefficients 


The   equation  actually  being  solved  is: 

AA 

AT  r   R 


(329) 


YAenf    analogous  to    (321): 

iX"^-      4       COAJD   (A/    1HL  (330) 

Ixl  |A\  ^   ; 

The  relative  error  in  the  residual  can  then  be  no  larger  than  the  rela- 
tive error  in  the  coefficients  in  A  times  the  cond(A). 

Finally,  an  approximate  accuracy  check  may  be  obtained  from  |j>J  : 

"  IE  or'kj  O  ,    , 

ij-    j^-j     ^  |0     ,  (331; 

then  X  is  probably  correct  to  p  digits. 

Equation  (300)  is  assumed  to  be  accurate  within  a  small  measurement 
error  and  equations  (303)  and  (30$)  are  considered  insufficient  to 
remove  the  inversion  instability.   There  are  several  techniques  avail- 
able to  remove  thses  instabilities  of  which  three  will  be  briefly 
covered. 

The  first  method  attempts  to  smooth  the  instabilities  by  one  of  two 

means.  Define: 
a(w)  =  (G-flM)T(G-AM)  +  yi(N-^)T(M-N0Vy2(6N)T^M)     (332) 


where   u(ll)  -  the  smoothing  measure  (a  number) 
}f   ,  $2  >  0,  numbers 

N«  -  a  priori  information  function,  vector 
B  -  matrix  describing  some  desirable  smoothing  of  II 
(G  -  AN)  (C— All)  -  a  measure  of  the  accuracy  with  which  N  satis- 
fies equation  (300) 

(IT  -  IL)  (N  -  N«)  -  a  measure  of  the  departure  of  II  from  the 
expected  II 

(BIT)  (BIT)  -  a  value  that  measures  the  departure  of  IT  from  ideal 
smoothness  (BN  =  0) 

The  smoothing  technique  involves  the  minimizing  of  u(ll).  Some  matrix 

differentiation  rules  will  be  needed  jj^j  : 

*X  -  /X  (333) 

~b%r     =  A  (334) 

<*XT  =  "1^  ax1"  (335) 

where  X,  Y,  Z,  are  vectors  and  A  is  a  matrix,  all  of  the  proper  dimensions 
such  that  the  indicated  products  are  defined. 

Condition  1:  o^    -   0;  requires  a  priori  knowledge,  N„. 
(332):    U,0V)=  yo00*  /,y,(M)  (336) 

where     V0(Aj)  -"  (G  -AA/rCG-AM)  (337) 

Vl(N)--(MX)T(Ni-Nia)  (33g) 


The  required  condition  is: 

U|(M)=  minimum    =>  ^q    U,M  1  =  0 


M-  M 


Utilizing  (333)'-.(  335 ) :    where   1    is  the  identity  matrix: 
^   =    -2/TCC-An) 


U 


(339) 

(340) 
(341) 


Use  (340)  and  (341)  in  (339): 

N,  =  (flTfl^,l)"'(flrG  +  y,/V0)  (3i2) 

This  is  the  required  solution  to  equation  (300)  given  the  constraint 
function  IL..   Note  that,  by  setting  5f,  =  0,  one  obtains  equation  (305). 
Since  the  multiplier,  V|  ,  represents  the  difference  magnitude  between  11 
and  lL,  one  might  expect  to  use  the  estimated  error,  (331 ),  as  a  first 
guess  at  J,  . 

Condition  2:  /,  =  0;  requires  no  pre-conceived  solution. 

(332):   u2(n)  =  y,(N)  *yz  y,M  (3«) 

where   yo  -  as  before  (337) 

y2f/v/)=  (bn)tCbmD  (3a) 

The  most  common  selection  for  B  requires  the  solution  to   be  smooth 
through  the  second  derivative.      In  difference  equation  form: 

Mj   =  bt(ri)  (345) 

tfj    =7?,-     Arj     =     7-^-^77     -    N0  (346) 

Assume  equally  spaced  quadrature  points: 

'I'^O-OtjES  (343) 

^rj  =  rjM"  rj    =     r*~r'    -    Ar  (349) 

Note  thatAr.  is  independent  of  ja   Thus: 

Ni  =aT*  (^-2A/j,,  +  Mifl)  (350) 

2 

Since  only  the  coefficients  of  K.  are  important,  ignore  the  l/^r  scale 

factor  and  write: 


(351) 


^  r 


or  N   -  BIT 

in  matrix  form  where  3  is  a  (k)x(k)  matrix  defined  ac 

/  1  -2     1      0  ...  0  0  0 

0  1    -2     1  ...  0  0  0 

0  0      1    -2  ...  0  0  Q 

0  0     0      1  ...  0  0  0 


(352) 


•      •      . 


•      •      • 


•      •      • 


•      •      . 


... 


(353) 


0  0  0...  1-2  1 
0  0  0  ...  0  1  -2 
0     0     0     ...      0     0      1 


The  vector  BIT  is  thus   the  discrete  analog  of  the  second  derivative  of  the 
unknown  vector  II.     Other  forms  for  3  may  be  used,   notably  reference  L20J  . 
As  before,    the   constraint  shall  be  to  minimize  the  quantity  u«(N)   in  an 
effort  to   achieve  the  smoothest  function  II  with  a  minimum  of  error  in  Gil: 

UtM*    minimum   =>    ?fi   ^MJm-A     =0  (354.) 

and  j^Yz  -     2BTBN  (355) 

Mow,    (340)   and   (355)   in   (354): 

n2  =  (a7\  +  ;t2btb)~V&  (356) 

This  is  the  required  solution  to  (300)  which  satisfies  the  constraints 
outlined  above.   Again,  setting  ^2  =  0  yields  (305).   A  recent  work 
into  this  particular  solution  has  yielded  considerations  for  the  choice 
of  ^2,.J20J    If  the  factor  ^  ^s  too  small,  the  instabilities  arc  not 
removed  from  the  final  solution.   And  if  chosen  too  large,  the  system 
becomes  overconstrained  and  independent  of  the  measurements.   The  best 
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choice  appears  tc   be    #2~10    '   and  should  be  so   chosen  until  further 
efforts  provide  a  better  choice. 

Both  solutions  via  this  technique  and  the  results  from  the  next 
method  may  be  improved  in  the  iterative  fashion  as  outlined  in  equations 
(321),    (325)- (327);    IIA.     The  algorithm: 

(initial  estimate)  1.     obtain  IV    '   by  (34-2),    (356),   or  (366) 

!   th         .,     t  ■v  0       0(m)       n        ,*(m)  r  T-n\ 

(m       residual)  2.     Rv       -  G  -  Mr    '  (357 J 

t   th  ..      ,    >  ~        ,-(m)  _   n(m)      £(m)      .-"!„(□)  ,,cm 

(m       error  estimate)  3.     Air    J  -  ?r    =>  ^v    '  =  A     R  (353) 

/      4..   j.  \  ,   *(m+l)  „  *(n)   ,  *(m)  m,m 

(new  estimate)  4..   ..      -  ^    +  NN  (35?) 

(mil 

(test)  5.  is  |R  'I  <  some  small  positive  number  $   ? 

no:  n  =  a  +  1,  go  to  step  2. 

A(m-H  ) 
yes:  done,  the  solution  is  N      to  an  accuracy  on  the  order 

of  $  . 

The  next  method  involves  the  computation  of  some  statistical  para- 
meters. Only  the  general  outline  shall  be  described.  See  reference  (_2lj 
for  details.  Let  "  and  /^   be  positive  definite  matrices  (all  eigen- 
values positive)  that  act  as  weighting  functions,  then: 

v^)--(G-ANyrxt>'Ati)+(.«-Ho)  +  r2w-Aj0)  (360) 

where  //  ^s  (n)x(n)  ^-d   /Q  is  (k)x(k). 

Let  Nn  be  the  mean  vector  for  II  as  compiled  from  previous  statistical 
data  (model).   Define  the  exoected  value  operator  E  as: 

E/jaooja  ^™li  <    3W+6O  (361) 

where  f(x)  is  the  probability  density  function  for  g(x).  Let  the  vector 
£  have  zero  mean  and  measure  the  error  in  G.  Define  the  covariance 
matrices  thus: 

SH-.E[(*-BHX«-e«r]  {362) 
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St     =    E&Cl  (363) 

Assuming  a  nornal  distribution  for  6  and  N,    then  the  Gauss-Markov 


A 


theorem  yields  a  statistically  optimum  estimate  IT     by  minimizing  v(ll) 


lor: 


P.  =  s" 


£  (364) 

T2    =     SJ  (365) 

then:           N3   =     X  ''A""  5^  (0  -  A  M0V  N0  (366) 

where           X     =      $'tl  +  A^ Sfc  A  (367) 

A  i 

The  resulting  covariance  matrix  for   (IT  -  II   )  is  found  to   be  X     .      Thus, 
the  expected  mean  square  error  may  be  written  as: 


e[(m-As)t(m-^]=  iracz  x 


(36S) 


Without  solving  equation  (300),  one  can  thus  obtain  an  accuracy  estimate 

A 

for  II  under  the  stated  statistical  conditions.  One  can  also  apnly  the 

A 

IIA,  (357)  -  (359),  to  enhance  the  estimate  U   . 

s 

The  last  technique  to  be  considered  is  an  iterative  scheme  somewhat 
analogous  to  the  iterative  improvement 

r  1  (o) 

the  Landweber  iteration. {j4J-   ■^n  SenersJ-»  sn  iterative  guess,  IT  '  is 

made  based  on  some  a  priori  knowledge  and  then  improved  successively 
to  form  a  sequence  whose  m  L    member,  I,   ,  transforms  to  the  measured 
G  as  m  tends  to  infinite.  The  techniques'  success  is  dv.e   to  the  auto- 
matic  constraining  of  the  oscillatory  nature  of  U         for  small  m.  This 
aliov/s  en  approximation  to  G  early  enough  to  avoid  the  oscillations  if 
the  sequence  is  terminated  when  the  residual  error  is  on  the  order  of  1 
error  in  G,  The  Landweber  algorithm  uses: 

NL  s  ML    +  AT(G-ANL    )  (369) 

where  js^j   reflects  the  expected  solution.  If  the  coefficient  matrix 
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-1       "U 


•las  aii.  positive  components,    tine  convergence  wi.^1   oe  improved   oy  de— 
fining  a  diagonal  matrix  D: 

D:  dj]  =  (|(a7   i-  I,...,k  (370) 

end  rewriting  the  estimate  as: 

N0  =Mo   +  t)ATCG-AM0    )  (371) 

Tliis  completes  this  short  survey  of  inversion  techniques.  A  common 
feature  lor  all  of  them  is  their  use  of  a  priori  knowledge  to  constrain 
the  selection  of  probable  solutions  in  order  to  enhance  the  choice  of 
the  proper  solution.   Indeed,  this  criterion  must  be  met  by  any  tech- 
nique the  experimenter  chooses  for  inversion,  if  he  is  to  have  any 
confidence  in  the  results. 


j^-i^O  -LVJ 


This  paper  has  endeavored  to  present  the  requirements  for  a  device 
designed  to  map  the  tropospheric  aerosol  properties  on  a  real-time  basis. 
The  specific  model  equations  were  derived  after  3.  theoretical  development 
was  performed  in  such  detail  that  the  user  became  familiar  with  the 
limitations  of  the  theory,   A  candidate  for  such  a  device  that  would 

r   chc 
..^  ijjj-iiii  system  allows  for  the  direct  selection  of 
the  hey  paramaters  of  the  theory  and  yields  as  output  the  required  in- 
puts for  the  data  inversion  to  begin.  A  few  inversion  techniques  were 
presented,  but,  due  to  the  large  available  extent  of  such  methods,  the 
subject  was  by  no  means  exhausted.  The  aerosol  IIP  A?,  mast  be  capable 
of  polarization  measurements  and  fast  wavelength  scanning  to  increase 
the  number  of  independent  observations  made.  Such  system  capability 
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would  decrease  the  required,  amount  of  a  priori  knowledge  the  experi- 
menter must  decide  on. 

Both  present  and  future  developments  in  fast  processors  and  versatile 
lasers  argue  strongly  for  the  design,  construction,  and  use  of  a  LIDAR 
aerosol  analysis  system. 


90 


APPENDED:   COMPUTATIONAL  FORMULAS  AIHD  ALGORITHM 

Tiais  section  presents  some  further  equations  to  facilitate  the  compu- 
tation of  the  required  fuctions  from  the  Mie  theory. 

Sachs catter:  ©  =  1C0° 

TTn  0M°;  =  (-I)""'1     $^  (372) 

Tn(/80°)  =.  -TTn  (l80°)  (373) 

5,  (190s;  ^  °|i  C-0n(n  +  l)(bn-An)  (374) 

5,(i809)   =   S,  Ciao0)   ■  (375) 

backscatter  gain:   &a  QstaO^O")  =^|  ^  C-0%+i)(t»n-a^  (376) 
Scattering;  functions: 

recurrence  relations: 

7n  (cos  e)  r  cos  q  TTn  (cos  e)  -  sin2efrn/Ccos  e)  (377) 

1T0(co59)=  CcosgX^t)^-,^  e)-(^)TTn,2Cco5e)  (378) 

tl/f^ose)  sf2n-0  fTn.,CcoS  ©;+  fl£2  (cos  ©)  (379) 

1T0  (cos  e)  =  O  (380) 

fh  (COS    ©)    =    I  (331) 

C(cos9)  =  0  (3S2) 

(COS  ©)  =    O  (3S3) 


%' 


Scattering  coefficients:  The  equations  for  a  and  b  may  be  written  in 
terms  of  the  logarithmic  decrement  function, ^ (Z) .  The  equations 
below  have  been  so  written  and  then  incorporated  into  an  algorithm  that 
was  used  to  compute  the  scattering  and  extinction  efficiencies  for  a 
particular  case  of  m  =  1 .  54-O-iO.COI )«  These  values  were  then  plotted 
on  the  enclosed  graphs  to  picture  the  oscillatory  behaviour  of  the 
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scattering  efficiency,  Q   ,  for  ck    from  0.1  to  17.0  in  0.1  increments. 

S  CO. 

The  data  points  were  connected  by  straight  line  segments  to  complete  the 
visualisation.   It  must  be  realized  that  further  structure  may  appear 
under  finer  resolution.  The  routine  was  programmed  for  the  Hewlett- 
Packard  model  4-1 G  machine  with  the  program  listing  accompanying  the 


.igures. 


(325) 
(386) 


Note  in  the  following  that  only  ( -.($/  requires  a  complex  argument. 


(103) 


(187) 


rnc^  rnSMhjr! 


yj  (<*.)  =  sin  ^ 
Xo  (<*)  -  cos  (<*) 

du)  %(<*)  =  y«w+ix„(^ 

^CV).       3n  CO  , 


v*c*o 


387 

383 

3S9 
390 
391 
392 
393 
394 

395 
396 
397 
393 
399 
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(400) 


/aCn-(31%,C/S)J  l„   W  ton/9 


Algorithm 


=>  Qw(-)  =  ^T  lif-) 
=>  Q«t(«)  =  £t%i.(<*) 

Begin:      A.      1.      Input:       n  0p£  ,  VC 
2.      compute:     mx    3    my 

B.  1.      input:  (S^cT    (required  accuracy) 

2.  compute:    /3*  ^  &  y 

3.  assign:     /?.-/  g,  C°0  =  g2  W  =  O 
4..      compute:        f   (<*)  ¥  (*) 

C.  I.      compute:      )"„' (".*. )  }    J£  „'  (<* ) 

2.  compute:       "n  W  0.o  («)t  bn  (*) 

3.  compute:        0^)  «  ?,  f-)  ♦  (2  „fl)   £/i„  Mj  '+  |  kn  (-  )  |  *} 

Cn)<;«)  =  g2  C«)  f  (z«*/ )  Re   {a.n<*)  t  bn  C«;  j 

yes:      go  to  D.I. 
no:        go  to  C.5. 
5.     test:     is       )  ^  (<*)  -  %  ,(*)  [  £    ? 
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2 

4-.      test:      is        i     (n) 


yes:      go   to  D.1. 
no:        go   to  E.I. 
D.      1.     assigr.:      <J(  O)  =    ^O) 

2.  compute:       /^f/   C^J  .,    ^n*/    '  °^ 

3.  assign:         ^_  ;   (<*)  =     ^   («*) 

n  =  n  +■  I 


z..     go   oO  G.  1 . 

GO 


/  2    \      tnj 

E.   1.   compute:   Qgca  £*)  =  («^y  ^ '  ^ 

2.   output:  n( number  of  steps) 

Q    («*) 
sea 

Qext  (or') 


Stop. 


94 


o 

"a, 
m 

U 

o 

Jh 

<B 

-P 
O 


0 

I — 

SI 

o 

•H 

o 

co 

• 

o 

CO 

•H 

— » 

I 

ta 

T— 

(4 

' * 

Q 

^f 

> 

U~\ 

• 

>> 

« — 

o 

c 

II 

o 

•H 

6 

O 

•H 

«H 

<H 

O 

M 

q 

•H 

k 

O 

-P 

.p 

cd 

o 

CO 

CD 

to 

•H 


tf 


M 
O 

-P 

0) 


o 

N 

•H 
00 


"SOS 

^    '^0U3T0TJJ8   Sirujao.q.tjoc 


95 


CD 
U 

o 

n, 

03 

U 

© 

-O 

b 


O 
N 

•H 

01 

CO 

3 

CI 
U 

o 

> 

o 
d 

© 

•H 
O 

(1) 

to 

•H 

?H 

0) 

-P 
■P 
cd 
o 

C/2 


O 

o 
o 

•H 
I 


cv 


0 
f-i 

•H 


m 

CD 

b 


o 

N 
•H 
00 


ft    'itouaxoTjja  Strr.iaq.q.120 


96 


;  ii  --  i 


0  1  ♦'_ SL   ' 

61  '  *HD£.S 

63  ri\r-~ 

66  PRO*' "  •" 

if.  I  j   I  K; 

,-.  n  n  .■•         -.  \ 

vO  P.  -  -     '.  - 

85  * 

18  C'rlc 

is  3 to  ;.-. 


!i* 


15  "HLF 

'-'- 

16  pROf 

I7  S'O 

I  1  \  '£ 

y-r 

19  'flCl 

lljRfiCY  ? 

28  PROMPT 

21  F3? 

5  - 

22  GTO 

91 

23  1  E- 

.  7 

24»LBL 

91 

25  STO 

95 

26  RCL 

64 

2?  RCL 

^2 

28  * 

il*  7  3>  »  U 

86 

76  RCL 

94 

31  RCL 

93  v 

32  * 

33  STO 

97 

34  8 

35  STO 

93 

36  STO 

V   1 
36  STQ 

89 

18 

84 

40  SI? 

41  STO 

i . 

43  COS 

4*  ST*.  1  & 

45  RCL  11 

46  RCL  64 
4?  f 

43  RCL  '3 


i1* 

*i7 


i"1   '.-y 


bo  biU 
64*LBL 
65  RCL 


7  8  1 


r  ; 

M 

- 

74 

STO 

t? 

1  -' 

1J 

76 

RCL 

14 

77 

RCL 

84 

!  0 

79 

DP' 

18 

88 

i 

31 

+ 

00 
0(. 

* 

83 

- 

34 

STO 

IS 

85 

RCL 

17 

86 

RCL 

12 

87 

/ 

83 

STO 

19 

39 

RCL 

18 

98 

RCL 

17 

91 

RCL 

14 

92 

RCL 

12 

93 

XEQ 

■  7  / " 

94 

STO 

26 

95 

XOI 

9b 

STO 

21 

...:'o  r,Ll 

ii£ 

161  RCL 

18 

ice  * 

103  RCL 

87 

184  RCL 

66 

185  XEQ 

"It-" 

166  RCL 

67 

107  RCL 

06 

168  XEQ 

•?A   L 

109  RCL 

'  6 

no  m 

111  - 

112  STO 

i£ 

113  XM 

114  STO 

33 

115  RCL 

■?7 

116  RCL 

"•o 

1  i  7  DPI 

Ci7 

1 ;  i  P'.L 

t1  * 

118  CHS 

11 '^  RCL 

86 

126  CHS 

121  XEQ 

"Z*a 

122  RCL 

16 

123  * 

124  RCL 

67 

125  RCL 

66 

126  XEQ 

■7*« 

4^7  v  r n 

itc  sty 

"1/2 

128  RCL 

33 

129  RCL 

75 

136  KEQ 

•Z*" 

131  STO 

24 

132  XO' 

i 

133  STO 

134  RCL 

63 

135  RCL 

19 

136  CHS 

137  * 

173  RCL 

62 

139  RCL 

19 

148  CHS 

141  * 

142  RCL 

o-. 

143  RCL 

24 

144  XEQ 

.z+. 

1  ■+  J   [».  L  L 

14 

<  j  •   riiV 

i  "; 

147  XEO 

£/ 

148  RCL 

12 

149  * 

150  STO 

J- 

15!  XO'i 

152  RCL 

i     ; 

153  * 

154  STO 

~  7 

<  cc  n,--i 

:  j.j  Kl-l 

21 

156  RCL 

Oj3 

157  RCL 

37 

158  CHS 

159  RCL 

62 

169  CHS 

161  XEQ 

■7*.u 

162  RCL 

25 

163  RCL 

24 

164  KcQ 

"Z+* 

165  XEQ 

■i/z 

166  RCL 

33 

167  RCL 

7° 

163  XEQ 

"Z*" 

169  STO 

26 

176  x<r 

171  STO 

97 

172  RCL 

25 

173  RCL 

24 

174  RCL 

87 

175  CHS 

176  RCL 

62 

177  CHS 

178  XEQ 

"2*" 

179  RCL 

19 

136  + 

131  RCL 

14 

132  RCL 

12 

183  XEQ 

"  z  ■■' " 

184  RCL 

12 

185  * 

136  STO 

70 

I  Of    :V  .  > 

t 

loo  fct/L 

12 

139  * 

196  STO 

11 

191  RCL 

25 

192  RCL 

24 

97 


,■>  1 


194 

!  95 

RCL 

82 

196 

CHS 

197 

X26 

-  _  *  '■ 

19'S 

RCL 

22 

tQq 

2* 

286 

X£ii 

*2+" 

2*?! 

XE6 

■  1  .•  7  - 

282 

RCL 

77 

267 

7'": 

26* 

XES 

"Z*y 

)|jC 

STO 

iC 

266 

X '. ) ' 

267 

STO 

>£ 

28  8 

X<  -'i 

289 

R-P 

216 

v* : 

2! ' 

370 

■j£ 

212 

RCL 

-7 

":  1  7 

rvr-i 

.5  .'_  w 

tb 

4  e-p 


21b 

Kl'L 

32 

217 

•*■ 

213 

RCl 

16 

219 

> 

228 

* 

OOj 

1 

t 

Oi7 

* 

224 

RCL 

68 

225 

+ 

226 

STO 

36 

227 

RCL 

00 

228 

RCL 

26  ' 

229 

+ 

238 

RCL 

16 

231 

0 

■'.'  7 0 

* 

7"-i 

C  •-'  ■-' 

1 

234 

+ 

235 

+ 

236 

RCL 

09 

237 

t 

233 

STO 

31 

239 

RCL 

36 

249 

RCL 

68 

241 

- 

•-1   A   --\ 

A  85 

243 

RCL 

2i 

244 

/ 

245 

RCL 

65 

246 

*<  =  T 

*) 

247  GTO  63 

24y  RCL  31 

249  RCL  ¥i 
.58  - 

251  ABS 

252  RCL  7! 

"'  c  7  f 

254  RCL  65 


^jb  ij! 


;T0  63 

257  GTO  64 
258»LBL  63 
259  RCL  36 
266  SIC  66 

261  RCL  31 

262  STO  69 

263  RCL  12 

264  RCL  18 

265  2 

266  * 

267  1 

268  t 

269  s 

276  RCL  64 
271  / 

Hi.   kll  li 

077  _ 

4.1  ■J 

274  STO  15 

275  RCL  12 
27b  STO  1! 

277  RCL  15 

278  STO  12 

279  RCl  14 
238  RCL  16 
col   c 

282  * 
233  1 
284  + 

235  * 

236  RCL  64 


'j  r-.n      f 

5.  0  1 

238  RCL 

13 

239  - 

298  STO 

16 

291  RCL 

14 

292  STO 

13 

293  RCL 

16 

294  STO 

14 

295  RCL 

25 

296  STO 

L.J 

297  RCL 

24 

298  STO 

OQQ  1 

'"' '"} 

C?7     I 

386  ST+ 

10 

361  GTO  22 

~~c.    !  ■■ 

382»L8L  64 

7  z  '     n    r. 

.'jo  r-p 

363  2 

757  RTH 

364  fCL  64 

358+LSL  'I* 

365  m 

YZ-*   R-F 

36-  / 

366  rDfi 

36"  RCL  36 

36!  RDH 

36:  + 

362  R-P 

3  2":<  STO  34 

v-'  7  '■■■■    V 

Tip  ? 

364  R DM 

31!  RCl  64 

365  * 

712  Y.12 

3c 6  RDH 

3l:  . 

367  - 

314  RCL  31 

368  Rt 

315  * 

369  P-R 

31b  STO  35 

776  RTH 

317  BEEP 

371+LBL  'It 

318  "STEPS  =  " 

772  XES  •!.' 

319  FIX  6 

373  XEQ  "2* 

328  CF  29 

374  RTH 

321  ARCL  16 

375*L8L  "TN 

322  AVIEH 

37o  2 

323  PSE 

377  + 

324  "CROSSECTIOliS- 

373  STO  36 

325  AVI  EM 

379  COS 

■icb  rot 

386  XOY 

327  RCL  35 

70  •  v 

323  RCl  34 

382  + 

329  SF  29 

383  STO  37 

336  FIX  7 

384  CHS 

731  "EXT=" 

335  ETX 

332  ARCL  Y 

386  RCL  37 

333  AVIEH 

337  EtX 

334  PSE 

333  + 

73c,  PSE 

339  2 

336  PSE 

398  •■" 

337  "SCA=- 

39!  + 

333  fiRCL  X 

392  STO  38 

339  AVI  EH 

393  RCL  37 

346  STOP 

394  EtX 

34!  GTO  66 

395  RCL  37 

342*LBL  "Z+" 

396  CHS 

343  XOY 

397  EtX 

344  RDN 

398  - 

345  + 

399  2 

346  REM 

466  <" 

347  + 

401  RCL  38 

348  Rt 

402  t 

349  RTH 

463  RCL  36 

358*LBL  "1/Z" 

464  SIH 

351  R-P 

485  RCL  22 

352  XOY 

406  / 

353  CHS 

467  RTH 

354  XOY 

463  END 

op 


LIST  OF  REFERENCES 


1.  Arfken,  G. ,  Mathematical  Methods  for  Physicists,  2d.  ed.,  Academic 
Press,  1970. 

2.  Gaputo,  B. ,  Leonard,  D.D.,  and  Guagliardo,  J.,  "LIDAR  for  Hire," 
Optical  Spectra,  U,  pp.  57  -  62,  1980. 

3.  Daily,  J.'./,  and  Harlenan,  D.R.F.,  Fluid  Dynamics,  2d.  ed. ,  Addison 
Wesley,  1973. 

4..   Deirmendjian,  D.,  Electromagnetic  Scattering  on  Spherical  Poly- 
dispersions,  Elsevier,  1969. 

5.  Gerald,  G.F.,  Allied  Mumerical  Analysis,  2d.  ed.,  Addison  Wesley, 
1978. 

6.  Grossman,  R.L.,  "The  Atmospheric  Environment,"  in  Remote  Sensing 
of  the  Troposphere,  ed.  by  7.E.  Derr,  university  of  Colorado,  1972. 

7.  Hall,  F.F.,  Jr.,  and  Ageno,  H.Y.,  "Absolute  Calibration  of  a  Laser 
System  for  Atmospheric  Probing",  Applied  Optics,  9,  pp.  1320  -  1824., 


1970. 


8.  Hecht,  E.  and  Zajac,  A.,  Optics,  3d.  ed.,  Addison  Wesley,  1976. 

9.  Hoel,  P.G.,  Introduction  to  Mathematical  Statistics,  4.th  ed.,  Wiley, 
1971. 

10.  I  shim  am,  A.,  Wave  Propagation  and  Scattering;  in  Random  Media,  2 
vols.,  Academic  Press,  1973. 

11.  Junge,  C.E.,  Air  Chemistry  and  Radioactivity,  Academic  Press,  1963. 

12.  Kerker,  M.,  The  Scattering  of  Light  and  Other  Electromagnetic 
Radiation,  Academic  Press,  1969. 

13.  Kolman,  3.,  Introductory  Linear  Algebra,  Macmillian,  1976. 

14-.  Landweber,  L.,  "An  Iteration  Formula  for  Fredholm  Integral  Equations 
of  the  First  Find,"  American  Journal  of  'MathonatJ.es,  73,  pp.  6 15- 
62/,,  1951. 

15.   Ott,  R.H.,  "Scattering,  Propagation,  and  Refraction  of  Electro- 
magnetic and  Acoustic  Waves,"  in  Remote  ''.'ens in ~  of  the  ^ropo sphere, 
ed.  7.E.  Derr,  University  of  Colorado,  1972. 


99 


1c.      Panofsky,    W.K.H.    and  Phillips,    II.,    Classical  Electricity  and 
Aarnetism,    2d.    ed.,    Addison  Wesley,    1962. 

17.  Schwars,    17. M,,    Intermediate  Electromagnetic  Theory,    Wiley,    1964-. 

18.  Schviesow,   R.L. ,    "Atomic,   Molecular,   Particulate,    and  Collective 
Generalised  Scattering,"   in  Remote  Sensing  of  the  Troposphere,    ed. 
by  *.-r.T.   Derr,   University  of  Colorado,    1972. 

19.  Colby,   S •!■!.,   Ctandard  ' Mathematical  Tables,    17th  ed.,    Chemical 
Rubber  Company,    1969. 

20.  The;;,    G.E.,    "Inversion  of  Optical  Scattering  and  Spectral  Extinction 
h'easurements  to  Recover  Aerosol  Size   Snectra,"      A"~S.ieA   Arties,    IS, 

p.   988  -  993,  1979. 


21. 


D    UX    CU  i  -.  ,         O   .  -  .  .         U-  -  a        uCOUniUSl)  -J.it., 

iluraerical  Solution  of  a  Tredholm  Aquation  of  the  Airst  hind," 

3    T    A    M       Traimfll     nf   'Ti irn pri  r»nl      lnalv<?iq       c      r>-n       2-7   -    2°^        1QAS 


22.  Strand,  O.IT.  and  Westwater,  E.R.,  "Inversion  Techniques,"  in  Remc 
ensing  of  the  Troposphere,  ed.  by  7. A.  Derr,  Tniversity  of 
do,  1972. 


juiui'ai 


23.  Strauch,  R.G.  and  Cohen,  A.,  "Atmospheric  Remote  Sensing  with  Laser 
Radar,"  in  Remote  Sensing  of  the  "Troposphere,  ed.  by  A.E.  Terr, 
University  of  Colorado,  1972. 

2A.   Van  de  Hulst,  H.G.,  Light  Scattering  by  Small  Particles,  Alley,  1957. 


100 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Information  Center 
Ganeron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  014-2 
Naval  Postgraduate  School 
Monterey,  California  93940 

3.  Department  Chairman,  Code  61 
Department  of  Physics  and  Chemistry 
Naval  Postgraduate  School 
Monterey,  California  9394-0 

4.  Professor  2.  C.  Crittenden,  Jr.,  Code  61  Ct 

(thesis  advisor) 
Department  of  Physics  and  Chemistry 
Naval  Postgraduate  School 
Monterey,  California  93940 

5o  LT  Dale  Robert  Ilamon,  DSN  (student) 
5S50  3.1/.  14th  Street 
Miami,  Florida  33144 


jooies 


101 


Thesis 

189505 

H1758 

Hamon 

c.l 

An  analysis  of  LIDAR 

atmospheric  remote 

sensing. 

HHT 


H      la 


■HE 
ISffi&K 


•■'.■. 


thesH1758 

An  analysis  of  LIDAR  atmospheric  remote 

III  III!  II  Hill 


3  2768  002  07600  2 

DUDLEY  KNOX  LIBRARY 


MP 

■"•'•'  V-X 

'■  -■ 

rS 

